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We propose a method that enables strong, coherent couphng between individual optical emitters 
and electromagnetic excitations in conducting nano-structures. The excitations are optical plasmons 
that can be localized to sub- wavelength dimensions. Under realistic conditions, the tight confinement 
causes optical emission to be almost entirely directed into the propagating plasmon modes via a 
\^ , mechanism analogous to cavity quantum electrodynamics. We first illustrate this result for the case 

' of a nanowire, before considering the optimized geometry of a nanotip. We describe an application 

f"**) ' of this technique involving efficient single-photon generation on demand, in which the plasmons 

, are efficiently out-coupled to a dielectric waveguide. Finally we analyze the effects of increased 

scattering due to surface roughness on these nano-structures. 
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I. INTRODUCTION 



In recent years there has been substantial interest in nanoscale optical devices based on local electric field enhance- 
ments and electromagnetic surface modes (surface plasmons) associated with sub- wavelength metallic systems. Surface 
" ^ ■ plasmons P| are electromagnetic excitations associated with charge density waves on the surface of a conducting ob- 
^-H [ ject. The unique properties of plasmons on nanoscale metallic systems have produced a number of dramatic observed 
P\| ' effects, such as single molecule detection with surface-enhanced Raman scattering (SERS) "5, 's'l, enhanced trans- 
04 mission through sub- wavelength apertures 0, , and enhanced photoluminescence from quantum wells 6] . There 
. is also considerable interest in these systems in applications such as biosensing , sub- wavelength imaging I3 ^ 
' and waveguiding and switching devices below the diffraction limit 0, Q ■ Such sub- wavelength waveguiding of 
plasmons in metallic nanowires has been observed in a number of recent experiments [TM [T3 . [l5| . 

At the same time, spurred in part by rapid developments in the fields of quantum computation and quantum 
information science, there has been strong interest in exploring new physical mechanisms that enable coherent coupling 
between individual quantum systems and photon fields. Such a mechanism would enable quantum information to 
+^ ■ be passed over long distances and long-range interactions between systems. These features are not only essential for 
^ . quantum communication [T^ Il7j but would also facilitate the scalability of quantum computers {l8[ . The required 
^ ' coupling between emitters and photons is difficult but has been achieved in a number of systems that reach the so-called 
Q^. "strong-coupling" regime of cavity quantum electrodynamics (QED) fT9lf2C}. l2l| . Recentlv several approaches to reach 



this regime on a chip at microwave frequencies have been suggested and experimentally observed (2ij |. 

which utilize coupling between emitters and modes of superconducting transmission lines. A key feature of these 
transmission lines is the reduction of the effective mode volume V^^ for the photons, which in turn results in a 



' substantial increase of the emitter- field coupling constant g oc l/^/V^g. Realization of analogous techniques with 
optical photons would open the door to many potential applications in quantum information science, and in addition 
lead to smaller mode volumes and faster interaction times. 

In this paper we describe a method that enables strong, coherent coupling between individual emitters and electro- 
magnetic excitations in conducting nano-structures on a chip at optical frequencies, via excitation of guided optical 
plasmons localized to nanoscale dimensions. The strong coupling occurs due to the sub-wavelength confinement and 
small mode volumes associated with the surface plasmon modes. We show that under realistic conditions optical 
emission can be almost entirely directed into these modes due to their large Purcell factors, in a manner analogous 
to cavity QED. We first examine the case of a cylindrical nanowire, a simple geometry where the relevant physics 
can be understood analytically, before considering the more optimized geometry of a conducting nanotip. We show 
that effective Purcell factors exceeding ~10'^ are possible in these systems, limited only by metal losses at optical 
frequencies. Because of these losses the plasmon modes themselves are not suitable as carriers of information over long 
distances. However, we show that the plasmon excitation can be efficiently out-coupled to a propagating photon by 
evanescently coupling to a nearby co-propagating dielectric waveguide, as illustrated schematically in Fig. ^ This can 
be used, e.g., to create an efficient single photon source, or as part of an architecture to perform controlled interactions 
between distant qubits. The achievable coupling between the plasmon and waveguide systems can be much stronger 
than the plasmon dissipation rates, and we find that single- photon generation efficiencies exceeding 95% are possible 
for the simple geometries considered here. 
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This paper is organized as follows. In Sec.^we calculate the mode structure of a conducting nanowire surrounded 
by some positive dielectric medium. We show that the nanowire supports one fundamental plasmon mode with 
significantly reduced phase velocity, which is tightly localized on a scale around the wire surface. We also 
calculate the dissipation rate of the fundamental mode as it propagates along the nanowire, due to metallic losses. In 
Sec, mil we calculate the emission properties of a dipole emitter near the nanowire as a function of emitter position and 
wire radius. We show that under certain circumstances, emission into the guided plasmon modes is greatly enhanced 
over decay into radiative and non-radiative channels. In fact, when optimized, the probability of emission into the 
plasmon mode approaches almost unity for small R and is limited only by dissipative loss of the metal. Because of 
its simple geometry, the nanowire is a system where the relevant physics can be understood and derived analytically, 
and from which we can proceed to design and understand better-optimized systems. In Sec. lIVI we consider one such 
system, a conducting nanotip. It will be seen that the enhancement of emission into plasmon modes found earlier is 
not exclusive to nanowires but arises quite generally as a feature of conducting nano-structures. However, we will show 
that the nanotip is an optimized geometry that can significantly reduce the effects of propagative losses even while 
preserving this enhancement. In Sec. we consider the problem of out-coupling the plasmon modes, and study in 
detail the interaction between the plasmon modes of our nano-structures and the guided modes of a nearby dielectric 
waveguide. We show that the plasmon modes can be efficiently out-coupled to the waveguide, and we propose an 
architecture for efficient single-photon generation on demand based on a tiered emitter/nano-structure/ waveguide 
system. We calculate the expected efficiencies for single photon generation, taking fully into account the propagative 
losses of the plasmons, the finite Purcell factors governing the interactions with the dipole emitter, and the non- 
unity coupling efficiency between the plasmon and waveguide modes. In Sec. IVII we consider the effects of possible 
imperfections to the system, in particular, the adverse effect of surface roughness on our nano-structures. In general, 
surface roughness can lead to radiative scattering of plasmons as well as increased non-radiative dissipation, which 
results in larger losses as the plasmons propagate along the structure. We calculate the effects of these two processes 
and find only moderate increases in the total loss under reasonable parameters. Finally, in Sec. IVIII we summarize 
our results, while outlining possible physical realizations and discussing possible future directions of research in this 
area. 



II. PLASMON MODES ON A NANOWIRE 



The method for calculating the electromagnetic modes of a nanowire is briefly outlined here, with details of the 
calculation given in Appendix ^ We consider a cylinder of radius R of dimensionless electric permittivity 62 , which 
is centered along the z-axis and surrounded by a second dielectric medium ei. While we are particularly interested 
in the case of a conducting nanowire surrounded by some lossless positive dielectric (Re 62 < 0, ei > 0), we note that 
at this point the discussion is quite general. Like any other simple geometry with a high degree of symme try, one 
can use separation of variables and find field solutions E, H to Maxwell's Equations in each dielectric region p5ll26j|. 
In cylindrical coordinates, the electric field is given by Ei(r) = £:i,„E,,„(fc,_L/^)e™'^e^'=ll^ where i = 1,2 denotes the 
regions outside and inside the cylinder, respectively. Here is the longitudinal component of the wavevector, which 
is related the vacuum wavevector fco = ^/c, electric permittivity ei, and transverse wavevector ki± by Cikg = fc^ + fcfx' 
and m is an integer characterizing the winding of the mode. A similar expression holds for the magnetic field H. For 
future reference, we also define the vacuum wavelength Aq — 27r/fco, and ki = y^fco as the wavevector in medium 
i. The coefficients Si^m and Tii^m multiplying the fields are not arbitrary but instead must satisfy a set of equations 
that enforces the necessary boundary conditions at the dielectric interface p = R. The existence of a non-trivial 
solution requires that the matrix corresponding to this linear system have zero determinant (detM = 0), which upon 
simplifying yields the mode equation |2^ , 

i?2 {kj^ kl^J \k2± J,n{k2±R) ki^ H„,{ki^R) J \k2± J^{k2±R) ki^ H^{ki^R) ) ' ^' 

One can use the above equation, for example, to determine the allowed values of k\\ as functions of m, i?, and e^. 

We now focus on the case of a sub-wavelength, conducting metal wire surrounded by a normal, positive dielectric. 
In Fig. 121 we plot the allowed wavevectors as determined through Eq. QJ, for such a system as a function of R 
for a few lowest-order modes in m. For concreteness, all numerical results presented in this paper are for a silver 
nanowire (or later, nanotip) at room temperature, Ag = 1 /im, and with a surrounding dielectric ei = 2, although the 
physical processes described are not specific to silver or to some narrow frequency range. The electric permittivity 
of silver at this frequency is e2~ — 50 -|- 0.6z, as given in |23|. In plotting Fig. 12 we have temporarily ignored the 
dissipative imaginary part of 62, although we will address its effect later. Ignoring Ime2 results in purely real values 
of k\\. 
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We first qualitatively discuss the important features of the plasmon modes illustrated in Fig. 13 before deriving them 
more carefully. It is clear from the figure that the longitudinal component of the wavevector exceeds the wavevector 

in uniform dielectric, fc|| > ki, which in turn causes the perpendicular component ki± = ^kf — /c^ = iKi± to be purely 

imaginary. Physically these relationships imply that the plasmon modes are non-radiative and are confined near the 
metal/dielectric interface, with the length scale of confinement determined by ^1/Kij_. Furthermore, these plasmon 
modes cannot couple directly to radiative fields, which have wavevectors k\\<ki. Of particular interest is the behavior 
of the plasmon modes in the nanowire limit In this limit, all higher-order modes (|m|>l) exhibit a cutoff 

as R^O, as derived in Appendix IbI while the m — fundamental plasmon mode exhibits a unique fc||cxl/i? behavior. 
Physically, in this limit, the rri — mode can be interpreted approximately as a quasi-static configuration of field 
and associated charge density wave on the wire. As such, R becomes the only relevant length scale, as the length 
scales l/|/ci| associated with electrodynamic behavior become unimportant. From the 1/R scaling of fc|| it follows that 
Ki±(xl/R, which states that the field outside the wire is tightly localized on a scale oci? around the metal surface. 
The corresponding small effective transverse mode area A^gocR^ is responsible for the strong interaction strength of 
the fundamental mode with nearby emitters, as will be discussed in following sections. We note that this behavior 
contrasts sharply with that of, e.g., a sub-wavelength normal dielectric waveguide or optical fiber, which runs into a 
"confinement problem" where the evanescent tails outside the device become exponentially large as R^O 0|. 

In practice £2 is not purely real but has a small imaginary part corresponding to metal losses (heating) at optical 
frequencies. Its effect is to add a small imaginary component to fc|| corresponding to dissipation as the plasmon 
propagates along the wire. In the inset of Fig. [3 we plot Re fc||/Im fc|| for the fundamental mode as a function of R. 
This quantity is proportional to the decay length in units of the plasmon wavelength Api=27r/Re fc|| . As R decreases, 
it can be seen that this ratio decreases monotonically but approaches a nonzero constant, as will be shown below. For 
silver at Aq = 1 yitm and room temperature and ei = 2 this constant is approximately 140. The fact that this ratio 
does not approach zero even as R—^0 is important for potential applications involving plasmons on nanowires, as it 
implies that the plasmons can still travel multiple Api for devices of any size. We also note that while all numbers and 
figures presented here are for room temperature, operating at lower temperatures might somewhat reduce the value 
of Im £2 due to decreased losses from phonon-assisted absorption . 

We now analyze the fundamental mode more carefully. For m — 0, one sees in Eq. ^ that one of the two terms 
on the right-hand side must equal zero. It can be shown that setting the first term to zero corresponds to a TE 
mode, while the other case corresponds to a TM mode (see Appendix 0. The TE mode equation does not have any 
solutions, and thus the fimdamental mode is a TM mode that satisfies the simplified equation [Toll30| 



k2± Jo{k2^R) kii_ Ho{ki±R) 
The fields themselves are given by (see Appendix ^ 



fci J^(fc2±i?) kl H'^ik^^R) _ ^ 



IK.\\ Z 



El = bii^-^^H;,{ki±p)p+^Ho{ki±p)zje 
E2 = 62('^^4(W)p+^Jo(W)5^"''"^ 



^2 



Hi = kii_hiH'^{ki^p)e' 

ujpo 



H2 = — fc2±62Jafe±p)e^'="^<^, (3) 

UJpo 



while the boundary conditions between the two dielectrics require that 

61 ^ k2± 4 {k2±R) 

62 ki^H[,{ki^_R)- 



(4) 



The 1/R dependence of k^\,ki± in the nanowire limit can be confirmed mathematically by considering Eq. in 
the non- retarded limit (c — > 00). In this case, ki± — ifc|| and the mode equation Q reduces to 

£2 ^ X^(fc||fl)/o(fc||fl) 

£1 Xo(fc||i?)/^(fc||i?)' ^ ' 

where Im.,Km are modified Bessel functions. The solution to Eq. jSjl requires fc||i? = C_i to be constant and proves 
the aforementioned scaling law for ku . It is also straightforward to see that when £2 acquires a small imaginary 
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component, the constant C_i becomes complex as well, and that Re fc||/Im fc|| takes on some fixed, non-zero value. 
No closed-form solution exists for the equation above, although when (corresponding to large |e2/ei|) the 

equation asymptotically approaches 



€1 (7-log2 + logC_i)(C_i)2' 

where 7«0.577 is Euler's constant. 

Finally, it should be noted that the components of in Eq. © are proportional to fcyfei^ or kfj^, while is 
proportional to ki±. Thus, in the nanowire limit when |fci^|cxl/i?, the magnetic fields are a factor of R smaller 
than the electric fields, which is consistent with this mode being roughly a quasi-static configuration. 



III. SPONTANEOUS EMISSION NEAR A METAL NANOWIRE 



The small mode volume associated with the fundamental plasmon mode of a nanowire offers ^a p ossible mechanism 
to achieve strong coupling with nearby optical emitters, in analogy to the methods of plll23.l23[ 24]. In this section we 
derive more rigorously the interaction between an emitter and nanowire, and show that under certain circumstances 
the small mode volume indeed leads to stron gly p referential spontaneous emission into the guided plasmon modes via 
a mechanism equivalent to the Purcell effect [sll in cavity QED. 

The spontaneous emission rate of a dipole emitter in general becomes altered from its free-space value in the presence 
of some dielectric body. In our system of interest, the dipole can possibly lose power radiatively to propagating photon 
modes, through excitation of the guided plasmon modes, or through non-radiative loss (heating) in the wire. The 
dipole in consideration can physically be formed by a single atom, a defect in a solid-state system, or any other system 
with a dipole-allowed transition. In Sec. IIII Al we calculate the radiative and non-radiative rates using a quasi-static 
approach, while waiting until Sec. IIII Bl to treat the plasmon decay rate more thoroughly. In Sec. IIII CI we show how 
the efficiency of emission into the plasmon modes can be optimized to yield Purcell factors in excess of ~5xl0^, and 
discuss the physical origins of this limit. 



A. Radiative and non-radiative decay rates 



In this subsection we derive formulas for the decay rates of a dipole near a metal nanowire into radiative and 
non-radiative channels. This calculation closely follows that of 32], but is briefly presented here for completeness. 

It is well-known that spontaneous emission rates can be obtained via classical calculations of the fields due to an 
oscillating dipole near the dielectric body [s^, and this method will be employed here. Specifically, we consider a 
(classical) oscillating dipole poe~*"* oriented along p and positioned a distance d from the center of the wire, and 
wish to calculate the total fields of the system. For nano-structures one can make a considerable simplification and 
consider the fields in the quasi-static limit (H«0) [sjl, which satisfy 

V-D = pext, (7) 
V X E = 0. (8) 

Here pext (r) is the external charge density. In the system of interest the external source is a dipole located at position 
r' outside the wire (with radial coordinate p' ~ d), which has a corresponding charge configuration 

Pext(r,r') = (po-V')<5(r-r'). (9) 

For simplicity we omit the harmonic time dependence from our expressions for the source and all fields. Note that the 
6{r — r') term above corresponds to a (unitless) point charge source, while the operator (po • V) generally converts 
the point charge solution to that of a dipole. It is therefore convenient to write E^ in similar form, 

E,(r,r')--V(po-V')<i>.(r,r'), (10) 

where 4>i(r, r') are "pseudopotentials" that satisfy V^<I>i = — (5(r — r')/eoei and V^<I>2 = 0. Here the indices 1,2 again 
denote the regions outside and inside the cylinder, respectively. Clearly, <I> physically corresponds to the potential 
due to a point charge at r', while the dipole potential follows from ^dip = (Po • V')<i>. 

To solve for the fields, it is convenient to further separate <I>i into "free" and "refiected" components 'I'o and ^r, 
respectively, where $r represents a source-free contribution that ensures that boundary conditions are satisfied and 
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$0 is the solution for a point charge in a medium of uniform electric permittivity ei. We will expand the known 
source term $o in basis appropriate for the cyhndrical geometry, and expand the source-free terms $^,2 in a similar 
basis that satisfies Laplace's Equation (V^$r,2 — 0). The unknown coefficients multiplying the basis functions of ^r,2 
will then be determined by enforcing the proper boundary conditions at the dielectric interface. Mathematically, the 
proper expansions are given by 



$o(r,r') - 



1 1 

47reoei |r - r'| 
1 

2^ 



5 y] (2 - <5m,o) COS {m{(f> - 0')) / dhcos {h{z - z'j) K,n{hp')I,n{hp) {p < p'), (11) 

$,(r, r') - I] (2 - <5m,o) cos (m(0 - 0')) / dh a,n{h) cos {h{z - z')) K^{hp')K„,{hp), (12) 

$2(r,r') = -_ y (2- 5,^,o) COS (m(0 -,/)')) / d/z /3™(ft) cos (/^(^ - z')) -^™(V)-^r„(V), (13) 



where am{h),f3m{h) thus far are unknown amplitude coefficients. We obtain a set of two coupled equations for 
am{h),Pm{h) by requiring continuity of $ and at the boundary, p = R. Because of the translational symmetry 
of the system, these equations are uncoupled in h and can easily be solved (this is in contrast to the case where 
translational symmetry is broken due to surface roughness, as discussed in Sec. IVIjl . The solutions are given by [33| 



{e-l)I'„^{hR)I„,{hR) 



eMhR)KiJhR) - €2Kra{hR)I'„^{hRy 



^ I„^{hR)K',^{hR) - K^{hR)r^{hR) 
' eiI^{hR)Ki^{hR) - e2K^{hR)lUhR)' ^ ^ 

where we have defined e = t^jti- Note that Eq. p4|l along with Eqs. (|ll|l - p3|) give the total electric field of the dipole 
and wire system. 

To calculate the radiative emission into free space, we consider the far-field properties of the system. Physically, 
the presence of the emitter induces some dipole moment (5p in the nanowire, which results in a total radiated power 
proportional to the square of the total dipole moment of the system, Pj-ad Trad |Po + (^Pp- We can determine 
(5p by finding the dipole-like contribution to the "refiected" potential ^dip,r = (Po • V')<i>r(r, r') far away from the 
source, which on physical grounds must behave like p~^ for large p. It is straightforward to show that the m = 1 
term in Eq. H12|) is responsible for this contribution, with all other m terms yielding faster decays in p. Because of 
the asymptotic behavior of Km{x)'^y/TT/2xe~^ when x^l, it can be seen that for large p the integrand appearing 

in H12|l is significant only over a small region h ^ p^^. As a result, we can safely replace Ki{hp') and ai{h) by their 
expansions around h — Q. After this simplification the integral can in fact be evaluated exactly and yields 



1 r 1 1 

r™"^' ~ -TT-o cos (0-0') / dhcoshiz- z')—Ki{hpf—h^R^ 

27r^eoei Jq hp' e + 1 



1 ^-1.^...^ P ^^^^ 

' V ( ra _L _ .n2^3/2' ^^^1 



47reoeie-Hl ^ ' p' (p2 + (^ z')2)3/2 ' 

with a corresponding reflected potential 

^ir.'^(r,r') = (po-V')$^=^Hr,r') 



■cos(0-0)— (16) 



47reoei e + 1 d'^ {p'^ + ^2)3/2 ■ 



In evaluating the above expression we have chosen p' — d and 2' = as the dipole coordinates, and po — PaP as the 

5p-r 

47reneir'^ ' 



dipole orientation. Comparing Eq. (|16|l to the potential due to a dipole Sp in uniform dielectric ei, V^p — — ^pjl 



we can readily identify 



<5P =Po— — r-Ts-P 17) 
e + 1 a 
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as the induced dipole moment in the wire, from which it follows that the radiative spontaneous emission rate is given 

by m 



*- rad 

To 



1 + 



1 <P 



{d>R) (18) 



Here Fq is defined to be the spontaneous emission rate in uniform dielectric ei (34] . Away from the plasmon reso- 
nance (ew — 1), the radiative decay rate changes slightly from Fq and reflects some moderate change in the radia- 
tive density of states in the vicinity of the nanowire. We note that this decay rate is well-behaved in either limit 
R, {d-R)^0. 

To calculate the other decay rates, one utilizes the fact that the total power loss of an oscillating dipole is proportional 
to the electric field in quadrature at the dipole's location, specifically, F^^q^^^j cxim (po • Ei(r', r')). Having divided up 
El into free and reflected components, we note that the contribution to Ei from the free field simply is associated 
with the decay rate in uniform dielectric ei, and thus we concentrate on the contribution from $r(r, r')- First we note 
that the coefficient ao{h) derived in Eq. p4|l contains a pole at the point where the denominator vanishes. This pole 
corresponds to an excitation of a natural mode (the fundamental plasmon mode) of the system. This can immediately 
be seen by comparing the denominator of ao to Eq. ((SJ, which gives the plasmon mode in the nanowire limit. The pole 
lies at hR — C_i and agrees with the plasmon wavevector derived in Sec.^ as expected. Evaluating the contribution 
of this pole to Er(r',r') gives the decay rate into the fundamental plasmon mode, and is discussed more carefully 
in the next subsection. At the same time, in the limit d^R one expects some type of divergence to occur in the 
non-radiative decay rate. Physically, such a divergence results from the large currents in the wire generated by the 
near-field of the dipole and their resulting dissipation. We can find the leading-order term to this divergent decay 
rate by carefully evaluating the leading-order divergence in the reflected field. 

In particular, while the pole associated with the m = term in yields the spontaneous emission rate into the 
plasmon modes, we will show that the mathematical origin of the divergence is the significant contribution to the field 
of an infinite number of terms with m > as d—^R. Specifically, in this limit, for a dipole oriented along p, 

rnon-rad _ 67reo Im p ■ E^(r', r') 
Fo 

^""^^ Imp- V(/5- V')3'r(r,r') 



o / dh K' (hd)^lm am(h) 

— ^ J2 / dhf„,ih,d,R). 



(19) 



The divergent nature of the above expression can be shown by examining the asymptotic behavior of the functions 



Uh,d,R)^< y + WV<l/ , (20) 



2d€i 



From the above expressions, we see that /,„ as a function of h has a characteristic width of about [2{d— R)] ^, 
yet at the same time the quantity m[R/d)^"^ reaches a maximum around to « 2(R-d) ^ d—>R. This confirms the 
non-vanishing contribution of an infinite number of terms with to > to the decay rate. The exact behavior of 
the functions fm at small h, including the peak around to, and the tails at large h is well-modelled by a Lorentzian 
approximation, 



Uh,d,R)^ -^- ,2.^ m2 ' (21) 
l + h^{d- R) 



which allows the integration and sum in Eq. H19() to be performed exactly. The resulting decay rate is given by 

Fnon-rad ^ , ^ ^ / £ 1 



l&klid-Rfel 
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Note that for |e|::^l and small Im e, Im (^f+r) ~2Im e/ (Re e)^, which makes it clear that the non-radiative spontaneous 
emission rate is proportional to the dissipative part of the electric permittivity. 

B. Decay rate into plasmon modes 

In this subsection we quantify the spontaneous emission rate Tpj of a dipole into the surface plasmon modes on a 
nanowire. As shown in the previous subsection, the coefhcient ao{h) characterizing the reflected field contains a pole 
at h = C-i/R that corresponds to excitation of the natural surface plasmon mode of the system. The contribution 
of this pole to the quantity Im(po • Ei(r',r')) yields the spontaneous emission rate into the plasmon modes and can 
readily be evaluated. Before proceeding further, we first note that in the presence of metal losses, the distinction 
between Tpj and Fuon-rad is not perfectly well-defined, since the plasmons eventually dissipate due to heating as well. 
Thus, for concreteness, we will define Tpj to be the decay rate resulting from the pole in the limit that Ime2 — 0, 
and take the plasmon wavevector fcy and C_i to be purely real in this subsection. In particular, for a dipole oriented 
along p, 

Tpi _ 67reo ^Imp ■ Er(r',r') 



67reo 



pole 

Im(p- V(p- V')$r(r,r')| 



= r4— Im( / dhh'^K^(hd)ao(h)] , (23) 

T^k^V^ \Jq J pole 

where we have explicitly indicated that we are interested in the pole contribution to the expressions above. It is 
convenient to explicitly separate out the pole of ao, and approximately describe the behavior around the pole's 
vicinity by 

,,,, _ 1 (ta-tijJilC-iUolC-i) 

"»<'•> :7(.-C_,/«)«^- 

where 

X{x) = eiIo{x)K'o{x) ~ eMx)I'^{x). (25) 

This separation allows us to easily evaluate Eq. ^Z',\\ and yields the decay rate 

KliC^^d/R) 
J^pl — apl^o- 



{koRf 

where we have identified Ki±~C-i/ R in the nanowire limit. The coefhcient apj is given by 

3(6i-e2)C^/i(C-i)/o(C_i) ^27) 



dx{C.,)/dx 



and most importantly depends only on ei^2- 

While the derivation above is straightforward, one can gain some physical understanding of the result and its relation 
to the Purcell effect by using Fermi's Golden Rule. This rule states that, once the plasmon modes are quantized, the 
decay rate is given by 

Tpi = 27rg\r,Lo)DiLo), (28) 

where g(r,Lu) is the position-dependent coupling strength between the quantized field and emitter, and D(lu) is the 
plasmon density of states on the nanowire. 

Canonical quantization of a dispersive medium is a difficult and subtle problem [s^ |3^ [sj IH, , and thus here 
we present a simple ad hoc quantization scheme that captures the relevant physics. To quantize the plasmon modes, 
we take the field solution in Eq. ^ and normalize the energy (again, ignoring Ime2) to 



d 



2 



hio^ I d''r[eo^iioe{r,iu)) Eir) + po H(r) ). (29) 
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The electric field term in Eq. H29|l gives the correct expression for the classical energy density in a dispersive 
medium [2^, and the coupling parameter at position r then simply follows through the relation g{r) — po • E(r)/?i. 
To evaluate the dispersive term, we assume that the conductor forming the wire exhibits Drude-like behavior with 
plasma frequency LOp, and that we operate well below the plasma frequency so that the permittivity is given by 
£2(1^) = 1 — Wp/w^w — iVp/ui"^. For such a metal this dispersive term is positive and given by ^ (^£2(0;)) «|e2(a;)|. 
Furthermore, we recall from Sec. ^ that in the nanowire limit the magnetic fields are smaller than the electric fields 
by a factor R. Combining these results, we find that the field energy is primarily electric, and 



huKieo / c;^r|£(r,w)| E(r) 



(30) 



Evaluating this equation readily leads to a normalization coefficient 61? 
sionless parameter that depends only on the permittivities £i_2, 



hukQeiR"^ /eoVCZiL, where F is a dimen- 
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(31) 



For a dipole oriented in the radial direction at position p' = d, the position-dependent coupling strength immediately 
follows, 
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Kiini^d). 



(32) 



The effective mode volume defined above is given by V^fj — VR^L and is proportional to the cross-sectional area of 
the wire and the quantization length L. This result reflects the transverse confinement of the plasmon on a scale 
comparable to R. Note that the presence of the l/VK;fT term in g is responsible for the strong coupling between 
plasmon modes and emitter as i?— >0. 

Assuming a Drude model, a scaling law for the density of states D{uj) — 2{L/2Ti){dk\\/dijj) can also be derived (the 
factor of 2 accounts for forward- and backward-propagating plasmons) : 



X Ldkn 

TT dUJ 
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(33) 



The important feature of Eq. (|33|l is the 1/R dependence due to the reduced group velocity doj/dk^^ocujR of plasmons 
on the nanowire. 

Combining the results of Eqs. H32f) and H33|l into Eq. H28|l . one finds that the decay rate into plasmons in the 
nanowire limit behaves like 



Fpi cx Fq 



KljK^^df 

{koRf 



(34) 



which agrees with the results derived previously. Again, the proportionality constant depends only on £1^2- Physically, 
the spontaneous emission rate into the plasmon modes increases like 1/R^ as i?— >0 due to the simultaneous reduction 
in group velocity {vgOcR) of these modes and an increase in coupling strength (g^cxl/i?^) due to the localization of 
the field energy to a region whose size is proportional to the cross-sectional area of the wire. 



C. Purcell factor of a nanowire 



Comparing the spontaneous emission rates given by Eqs. (|18|l . (|22ll . and (|26|l . we now qualitatively discuss the 
behavior one should expect as the position of the emitter is varied. In the limit that d/i?3>l clearly the spontaneous 
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emission rate is dominated by radiative decay and is equal to the spontaneous emission rate Tq in a uniform dielectric. 
As one brings the emitter closer to the wire surface, the change in the electromagnetic mode structure near the wire 
results in some modified radiative decay rate Fj-ad which never exceeds approximately 4ro for large |e|. When the 
emitter position d approaches (i~l/|fcij^|~i?/|C_i|, the emitter starts to be able to decay into the localized plasmon 
fields, with the rate scaling with wire size like 1/ . The spontaneous emission rate into plasmons continues to grow as 
the emitter is brought even closer to the wire edge, d^R. However, the efhciency or probability of plasmon excitation 
eventually decreases due to the large non-radiative decay rate experienced by the dipole very near the wire, which 
diverges like l/{d — R)^. We thus expect some optimal efficiency of spontaneous emission into the plasmon modes to 
occur when the emitter is positioned at a distance 0{R) away from the wire edge, and for this optimal efficiency to 
improve as R^O. 

This result is illustrated in Fig. where we have numerically evaluated the spontaneous emission rates derived 
previously. Specifically, we plot as a function of R the "error" probability Pe = 1 — rpi/(rpi + F') that a single, 
excited quantum emitter fails to decay into the plasmon mode. Here, F' = Fj-ad + Tnon-rad denotes the total emission 
rate into channels other than the fundamental plasmon mode, and the error probability has been optimized over the 
emitter position d. It can be seen that as i?— s-0, the probability of emission into the plasmons approaches almost 
unity. Examining this limit more carefully, the error in fact approaches a small factor Pgcxlm e/(Re e)^, explicitly 
indicating that the efficiency is limited by dissipative losses, as will be more carefully shown below. For the chosen 
parameters the probability of emission into the plasmons is well over 99% as i?— >0, with a corresponding effective 
Purcell factor Fpj/F' « 5.2 x 10^. Again, we emphasize that these properties are specifically a result of the conducting 
properties of the nanowire. This can be contrasted with emission into the guided modes of a sub-wavelength optical 
fiber, which drops exponentially as R^O due to the weak confinement of these guided modes (s^. In Fig. we 
plot logj^Q Pe as functions of R and d/R. It can be seen that achieving a large Purcell factor does not depend too 
sensitively on the emitter position d. 

We now prove that the maximum efficiency as i?— >0 is indeed limited by a small factor related to the dissipative 
losses of the metal. We consider the quantity 



T' (1 + RVd^Y + anon-rad(fco(rf - R))-' 

Fpi api{koR)~^Ki{Ki±d)^ 



(35) 



where ttnon-rad ~ 

(3/8ei/^)Ime/(Ree)2 is a small parameter explicitly characterizing the losses in the metal. Defining 
y={d — R)/R, and using ki±_k.C-i/ R along with the asymptotic expression Ki{x) « 7r/2a;e~^ for large x, we can 
re- write F'/Fpj as 



Fpl TTttpl 



+ 
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{i + yf 



Q^non-rad 
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(36) 



Note that the first term in the brackets corresponds to radiative decay and vanishes in the limit that i?— s-0, and 
thus the ultimate limit to the efficiency of plasmon generation is due to a balance between the plasmon and non- 
radiative decay rates. In this limit, a straightforward calculation yields a minimum in the expression above at 
yo = {1 — C-i + \J\ + C-i(4 + C_i))/2C_i, which confirms that the optimum position of the emitter is on the order 
of a few radii away from the wire edge, while the corresponding value of the minimum is proportional to ttnon-rad- 



IV. SPONTANEOUS EMISSION NEAR A NANOTIP 



In SectionsHTland llTll we derived and discussed the physics of plasmon modes on a nanowire and spontaneous emission 
of a nearby dipole emitter. For this simple geometry it was possible to find analytical solutions and understand 
the relevant physics of emitter/plasmon coupling in conducting nano-structures. In particular, it was seen that 
for such structures, the tight transverse confinement of the plasmon modes leads to a large effective Purcell factor 
for an optimally positioned dipole emitter as the relevant size scale decreases, with the maximum enhancement 
limited by non-radiative decay. At the same time, however, it is evident that the i?— >0 limit is accompanied by 
enhanced losses as the plasmon propagates, due to the tighter confinement of fields in the metal, and a reduction 
in the plasmon wavelength Api that could make out-coupling more difficult. Such factors could clearly impose limits 
for applications such as quantum information, but can be circumvented with simple design improvements. In this 
section, we investigate one specific design, a metallic nanotip. As in the nanowire case, one expects a sub-wavelength 
plasmon mode volume, determined here by the tip curvature, and an associated enhancement of emission into the 
plasmon modes. At the same time, though, the tip can rapidly expand to larger sizes where the propagative losses 
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of the plasmons are less severe, and where Api is not as smalL In the nanotip case we are not able to obtain full 
electrodynamic solutions for the plasmon modes. However, in a manner similar to that described in Sec. Illll we will 
calculate all of the relevant decay rates in the quasistatic limit and describe an approximative method to calculate 
the effects of propagative losses along the nanotip. We will also compare these results to those obtained via fully 
electrodynamic numerical simulations, and we find that these two approaches agree closely. 

In the following we will consider a nanotip whose surface can be parameterized as a paraboloid of revolution with 
symmetry along the z-axis. Specifically, we suppose that the surface of the nanotip is described by 

U-' + y' ^..\ (37) 



a paraboloid of revolution with apex at z = — Uq/2 (the reason for the offset of the apex will become apparent below). 
We now introduce a change of coordinates, 

X = uvcosfj), (38) 
y = uvsmcj), (39) 

z ^liu'- v') . (40) 

While these coordinates may seem awkward (note, for example, that u, v have units of ^length), they are convenient 
for deriving expressions for the fields and spontaneous emission rates, which we will then express in more "natural" 
coordinates at the end of the calculation. In these parabolic coordinates, the nanotip profile of Eq. H37|) is defined by 
a surface of constant v = vq. More generally, any constant v defines some paraboloid of revolution in this system, 
while the unit vectors u and v run normally and tangentially to these surfaces, respectively. 

Now, as in the nanowire case, we are interested in seeking the quasistatic field solution for a point charge source in 
the vicinity of the nanotip, from which we can obtain the field due to a dipole po. In particular, we seek solutions of 
the total field of the form (|10|l with appropriate boundary conditions. Like before, we separate the pseudopotential 
<&! outside the nanotip into its free and reflected components $o,r, and use an integral representation of the free 
pseudopotential suitable for parabolic coordinates, 

$o(r,r') = ^ ^ (2 - (S„,o) cos micj) - cj)') / dq qJ„,{qu)J„,{qu')I^{qv)K,n{qv'). {v < v') (41) 



27reoei 

m— 







Because $o fully accounts for the point source, $^,2 then satisfy Laplace's Equation. Using separation of vari- 
ables, it is straightforward to show that the solutions to Laplace's Equation are given in parabolic coordinates by 
^Jm{qu)Gi^m{qv)e^™''^, where Gi,m(gu) = Km{qv) and G2,m{qv) = Im{qv) are non-divergent functions in their regions 
of applicability. We then define the following expansions. 



oc 



$r(r, r') — (2 - (5„i,o) cos to(0 - 0') / dq qa„i{q)Jjn{qu')K„i{qv')J,n{qu)K,n{qv), (42) 

STreoei Jo 



OO 



$2(r, r') = — 2^ (2 - dm,o) cos m(0 - 0') / dq q(3,n{q)J,n{qu')K^{qv')J„,{qu)I„,[qv), (43) 

m— 







where the coefficients a,/3 will be determined by imposing boundary conditions at the nanotip surface v — vq. It can 
be easily shown that the continuity of $ and Dx imply that 

^ ^ (ei -e2)/;,(g«o)J,„(gt;o) ^^^^ 

(^2lmiQVo)K,n{qvQ) - eiI„i(qvo)K'^{qvo) ' 
^ . ) ^ -£i iImiqvo)K'^iqvo) ~ I'm{qvo)K„,{qvo)) 

Note that the coefficients am{q), along with Eq. H42I) . completely determine the reflected field. 

The calculation of the radiative and non-radiative spontaneous emission rates proceeds in the same manner as the 
nanowire case. To calculate Fj-ad, we again look for a dipole term in the far-field (large v) that corresponds to an 
induced dipole moment (5p in the nanotip, and then use the relationship Fj-g^doclpo + <5pp. At the same time, we look 
for a divergent contribution to the reflected field at the dipole location as its position v' approaches wq, which yields 
the leading term of the non-radiative decay rate through rnon_j.adC<Ini(Po • E,.(r',r')). This divergence is physically 
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due to the dissipation of divergent currents induced in the metal by the dipole. For a dipole positioned along the 
z-axis {u' = 0), 



Po-E,(r',r') 
Po •E,(r',r') 



27reoei 



dq -pTai{q)Kf{qv'), (po -L 
q^ 

dq -p^ao{q)Kl{qv'), (po || z) 



(46) 



as shown more carefully in Appendix El Mathematically, the divergence as v'^vq occurs due to the presence of a 
long tail ~e~^''^" """^ in the integrand for large q. Because of the similarity of the decay rate calculations with those 
in Sec. IIIII we simply state the results here, while providing more details in Appendix |0 For a dipole positioned 
along the z-axis at v ~ v' , the radiative and non-radiative spontaneous emission rates are given by 
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(47) 



and 
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(48) 



Finally we consider the decay rate into the fundamental plasmon mode of the nanotip, which is associated with 
the contribution of the poles in the integrand of Eq. (|46l) to Im (po • Er(r', r')). The presence of a pole indicates the 
excitation of a natural mode of the system. Examining the solutions to ao,i given in Eq. H44|l . one finds that ai 
has no pole in the range < q < oo. Physically, the absence of a pole means that a dipole simultaneously oriented 
perpendicular to z and located along the z-axis does not excite the fundamental plasmon mode of the nanotip. This 
is easily understood since a dipole oriented this way is anti-symmetric with respect to 180° rotations about z, while 
the plasmon mode is symmetric. On the other hand, ao does have a pole corresponding to plasmon excitation. This 
pole is located at qo — C-i/vq, where C_i is the solution to Eq. ((S)). Evaluating the contribution of this pole to the 
field is straightforward and yields a plasmon decay rate 
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e2)/i(C_i)/o(C_i) 
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(49) 



where x(a;) is defined in Eq. H25() . As in the nanowire case, the decay rate Fpi into the plasmon mode given by 
Eqs. (|49|1 and H25|) is evaluated in the limit that Ime2 = 0, such that go and C_i are purely real. 

Having derived the decay rates in parabolic coordinates, we now define a more natural set of parameters to describe 
the system. Let us introduce a length scale w that characterizes the nanotip via p{z) = \Jwz (z > 0), where p is the 
radius of the nanotip at position z (note also the corresponding shift in the apex of the tip from z = — Uq/2 to z = 0). 
Furthermore, let z = —d < be the position of the emitter [d is the distance between the emitter and end of the 
nanotip). In terms of these parameters, the spontaneous emission rates derived above can be re- written as 
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where a.^\ only depends on t\^2 and is given by 
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Having obtained the spontaneous emission rates near a nanotip into the different possible channels, it is once again 
possible to optimize the efficiency of emission into the plasmon modes for given w by varying the emitter position d. 
The corresponding optimized error probability, Pe = 1 — rpi/(rpi + T'); plotted as a function of w in Fig.|Sli. In 
analogy to the nanowire system, a large effective Purcell enhancement of rpi/r'«2.5 x 10'^ arises as >0, due to a 
balance between the small mode volumes associated with the plasmons and the non-radiative decay rate. In Fig. Et, 
we plot logj^Q Pe as functions of w and d/w. Once again, it can be seen that the error is not too sensitive to the 
emitter position. 

Because the plasmon modes here were obtained through a quasistatic approximation, this calculation yields no 
information about dissipative losses as the plasmon propagates along the nanotip. For example, in this limit H«0 
so one cannot obtain the Poynting vector for the system. To estimate the effect of propagative losses, however, we 
can make an eikonal approximation 40], assuming that the plasmons are emitted completely into the end of the 
tip {z — 0), and that the propagative losses thereafter at any position z are described locally by the nanowire solution 
at radius p{z). This motivates us to define an effective decay rate 

/ rziR) \ 

rpi(i?) = Tpiexp ( -2 Im fc,, {p{z))dz \ , (54) 

which equals the rate of decay into the plasmons multiplied by the probability that an emitted plasmon will propagate 
without dissipation to some final nanotip radius R. Here Fpi is the decay rate for the nanotip obtained earlier, while 
fc|| (p) is the plasmon wavevector for a nanowire of radius p. One can also define a corresponding effective error 
probability Pe{R) = 1 — rpi(i?)/(rpi + F') for the nanotip. Physically, Pe is the probability that the plasmon mode 
is either not excited by the emitter, or is excited but fails to successfully propagate to final radius R. In Fig. we 
plot this quantity as a function of R, when optimized over possible tip parameters w and d. It can be seen that the 
effective error probability for a nanotip compares favorably to that of a nanowire when koR ^ 0.05. In other words, 
the nanotip configuration is able to simultaneously exhibit a strong Purcell effect and reduce propagative losses. We 
note that the nanotip system also has the added benefit of generating guided plasmons along a single direction of 
propagation. 

Finally we discuss the limits of validity of the equations derived above for the nanotip. The quasistatic decay 
rates are valid in the regime |A:i|c?^l, which implies that propagative phases associated with the electrodynamic 

solution can be ignored over the length scales of interest. At the same time, Eq. (|54f) assumes that the plasmon 
mode at some tip radius p adiabatically follows the nanowire solution of corresponding radius. One can define 
an adiabatic parameter /3 = d(l/Re fc|| (/9))/c?z associated with the propagation, which must be small for such an 
assumption to be valid. Physically, a small P corresponds to a small correction to the propagative phase due to the 
variation of Re fc|| {p{z)), compared to the ~27r phase acquired over a distance of the plasmon wavelength. Assuming, 
for example, that we are considering sufficiently small length scales that Re fc||«|/c|| the 
condition that /3<d implies that the eikonal approximation is valid only in regions where z3>if/|C^i|. It can be 
seen that Zc=w/\C'^i\ represents some cross-over value, below which Eq. (154(1 clearly is invalid. On the other hand, 
ll/^ol = ^o/|C'^iH^/l^-il ssts the relevant length scale for the plasmons on the nanotip, and one expects dissipation 
to occur on length scales much longer than this. Thus, as long as the losses predicted by Eq. 154|) for z < Zc remain 
small, one can effectively use this equation for all z even if it is not strictly valid for z < Zc- A straightforward 
calculation confirms that the predicted loss, 1 — Fpi(i?c)/rpi^l — exp(— ImC_i/|C-i|), is indeed negligible, where 
Rc=R{zc). Finally, in practice, for the applications of interest we will be primarily interested in nanotip devices whose 
radii do not grow indefinitely, but rather expand until they reach some final, constant radius R. Such devices, for 
example, are more likely to be easily out-coupled, as discussed further in the next section. For Eq. ((54|l to remain valid 
for such a system, the spontaneous emission rate into plasmons for this device must be close to the rate calculated 
for an infinite, perfectly paraboloidal tip. This imposes the additional requirement that the final radius R be much 
larger than i?c- 

To check the analytical results derived above for the nanotip, we have also performed detailed numerical simula- 
tions using boundary element method (BEM) Details of our implementation are given in Appendix IdI BEM 
simulations are fully electrodynamic solvers of Maxwell's Equations, and they were used to obtain the classical elec- 
tromagnetic field solutions of an oscillating dipole emitter poe^*'^* near a nanotip. The results of a few sample 
simulations are shown in Fig. 0] for a tip curvature parameter kow = 0.022, final radius koR — 0.3, and varying 
emitter positions kod — 0.002,0.2,0.7. In Fig. we plot the quantity |Re(E x H*)|, which is proportional to the 
Poynting vector and corresponds to the total energy flux of the system. The nanotip is assumed to be composed of 
silver in a surrounding dielectric ei — 2, and its boundary is given by the dotted line, while the emitter positions are 
denoted by the circles. The total spontaneous emission rate is given via Ftotal = (^pi -I- r')cxlm (po ■ Ei(r',r')) and is 
determined numerically for each configuration by finding the total field at the dipole location. On the other hand, the 
effective emission rate rpi(i?) into the plasmons is determined by a best fit of the fields in the region of constant R 
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to the known plasmon solution on a nanowire given in Eq. (O, and then calculating the total power transport of this 
best- fit mode through the integrated Poynting vector. This total power is directly proportional to rpi(i?). The figure 
confirms the qualitative behavior that we expect and have described previously. In particular, the generated plasmon 
field and total spontaneous emission rate are largest for very small separations and decrease as the emitter is placed 
further away from the end of the nanotip. In Fig.^jD, we plot |Re(E x H*)|/rtotal; which is proportional to the energy 
flux normalized by the total power output of the emitter. This quantity yields information about the efficiency of 
decay into the various channels. For small separations {kod — 0.002), the plot is mostly dark, which indicates that the 
decay of the dipole is predominantly non-radiative. For kod = 0.2, the maximum (corresponding to bright spots in 
the plot) is located along the entire surface of the nanotip, which indicates highly efficient plasmon excitation. Here, 
although the total emission rate into plasmons decreases from the kgd = 0.002 case (as seen in Fig.0^), the efficiency 
increases dramatically due to less competition from non-radiative decay. Finally, for kod = 0.7, the maximum appears 
as the typical lobe pattern associated with radiative decay. 

In Fig. 13^, we have plotted the numerically optimized values of Pe{R) = 1 — rpl(-R)/(rpi -I- F') for a few values 
of R. It can be seen that the values of Pe obtained through analytical approximations and numerical BEM closely 
agree. Unlike the theoretical predictions, however, the numerically calculated error probability does not increase 
monotonically with R. We believe that the origin of this is that for the numerically optimized parameters, the 
condition i?3>i?c under which the theoretical predictions hold is only weakly satisfied, and the excitation region for 
the plasmons cannot strictly be thought of as a single point at the end of the tip (z = 0). 

V. SINGLE PHOTON GENERATION VIA COUPLING TO DIELECTRIC WAVEGUIDE 

We have shown in previous sections that a single emitter can spontaneously emit into the guided plasmon modes of 
a nearby nano-structure with high probability. This prospect of efficient conversion between excitation of the emitter 
and a single photon has a number of applications in the fields of quantum computing and quantum information. In this 
section, we consider one particular application, involving the use of such a system as an efficient single-photon source. 
The concepts behind single-photon generation on demand with an individual emitter in a cavity have been discussed 
elsewhere 42ii 4^ 44] and will not be presented in detail here. We note also that the ideas behind single-photon 
sources can be extended to create long-distance entanglement between emitters, as detailed, e.g., in jS^. 

Because of dissipative losses in metals, the plasmon modes are not directly suitable as carriers of information 
over long distances. We show, however, that plasmonic devices can serve as an effective intermediate step, and in 
particular can be efficiently out-coupled to the modes of a co-propagating dielectric waveguide. The single photon 
device is illustrated schematically in Fig. ^ In Fig. an optically addressable emitter with multiple internal levels 
sits in the vicinity of a conducting nanowire. The emitter is strongly coupled to the nanowire, such that single photons 
on demand can be generated with high efficiency in the plasmon modes by external manipulation of the emitter. The 
addressability of the emitter along with the internal levels allows for shaping of this single-photon pulse as 
illustrated in Fig. QJd. Here, a three-level emitter is shown with two ground or metastable states |s), which both 
have dipole- allowed transitions to the excited state |e). We assume that the system is prepared initially in the state 
|s), and that the |e) — \g) transition is coupled via the plasmon modes of the nanowire, i.e., the state |e) can decay at a 
rate Fpi into \g) by emitting a photon into the plasmon modes. In addition, there is a small rate F' at which the excited 
state can decay without emitting a plasmon. A single photon in the plasmon modes of the nanowire is generated 
with high probability by exciting the \s) — \e) transition with some external pulse Q,{t) and the subsequent decay into 
\g), with the shape of the single photon wavepacket controlled by the shape of il(t). We further assume that the 
plasmon is then evanescently coupled to a nearby dielectric waveguide, as shown in Fig.^, which co-propagates with 
the nanowire over some distance L^x over which this coupling is non-negligible. The coupling is a reversible process, 
and the distance L^x is optimized to maximize efficiency of ending up with a single photon in the waveguide (i.e., 
to prevent further Rabi oscillations back into the nanowire). A similar setup with a nanotip is illustrated in Fig.^. 
Here the nanotip radius p{z) expands to some final radius R at which point coupling with the waveguide starts to 
occur. Initiating the coupling once the nanotip has reached a constant radius allows the two systems to be easily 
coupled, as discussed below. When optimized, we estimate that single-photon generation efficiencies exceeding ^--^90% 
are possible in this tiered configuration. 

To treat the problem analytically, we consider the simple situation of our nano-structure coupled to a cylindrical 
dielectric waveguide {e.g., an optical fiber) of radius Rg, such that the modes can be calculated analytically using the 
methods described in Appendix ^ It can be shown that the fundamental modes of the waveguide are degenerate 
m = ±1 modes that are not cut off as Rg^O. The dependence of their wavevector fcy on Rg is shown in Fig. [SJ 
for a core permittivity = 13 and surrounding permittivity ei — 2. These parameters correspond closely to that 
of a Si/Si02 guide at Aq = 1 /xm. To simplify the calculation, we also assume that coupling between the wire and 
higher-order waveguide modes is negligible. This can be achieved, for example, by operating below the cutoff radius of 
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higher-order modes or by operating with sufRciently large wavevector mismatch between the plasmon and higher-order 
guide modes. 

We make the ansatz that the total field of the system is given by a superposition of the unperturbed modes of the 
nano-structure and waveguide. While this cannot strictly be correct, as such a solution violates boundary conditions 
at each interface, we rely on such an assumption to give us the correct qualitative behavior without resorting to more 
complex numerical calculations. Specifically, we assume that the total electric field for the system takes the form 

where /i indexes the nano-structure {w) and waveguide {g) systems, and i = 1, • • •, iV^ runs over the modes of system /i. 
In the following we will explicitly treat the nanotip case, where the plasmons propagate in a single direction, although 
this argument can easily be extended to the nanowire. We emphasize that we are considering coupling of the plasmon 
mode to the waveguide once the nanotip has already expanded to its final radius i?, at which point the plasmon mode 
solution becomes identical to that of a nanowire. In our case, = 1 as we only consider the fundamental plasmon 
mode of the nanotip, while Ng = 2 as we take into account the two degenerate, co-propagating fundamental modes 
of the waveguide. E^_i(r) here represents the unperturbed solution of mode i in system /i (without the presence of 
the other system). A similar expression holds for the total magnetic field. 

With the ansatz of Eq. (|55ll for the total field of the combined waveguide and nanotip system, one can derive 
exact equations of evolution 47] based on Lorcntz reciprocity for the coefficients C^,i. Explicitly separating out the 
plane-wave dependence of the unperturbed fields, ^(r) — E^ ^(p)e*'^iif the -t- Ng coupled- mode equations 
take the form 

E E^M.;..^w^ = *^^o E E^M.....wa,,w, (56) 

u—w,g j — 1 u—w.g j — 1 

as derived in detail in Appendix IeI The coefficients to the system of equations above are given by 

P^,.;.,,(^) = e^('=ii--'=ii;.)^y' dp (E,,,(p) X H;,,(p) +E;,,(p) X H,,,(p)) . z, (57) 

if^,,;.,,(^) = e'('=ii--'=ii:.)^y' dpE,,,(p).E;^,,(p)(eT(p)-e.(p)), (58) 

where erip) is the electric permittivity of the combined system. Clearly, the presence of the phase factors 
i)2 equations above indicate that, at least under weak coupling, significant power transfer be- 

tween the two systems will not take place unless the two systems are approximately mode-matched with respect to A;|| . 
In practice, this implies that for a final tip radius R, there is some ideal waveguide size Rg that allows for maximum 
transfer efficiency between the two systems. A similar optimization of the waveguide parameters exists in the case 
of arbitrary coupling strength between the two systems, although this problem is more complex because one must 
account for factors such as the phase shift of one system due to the other. We emphasize that the coupled-mode 
equations above are exact within the ansatz of Eq. H55|) . For example, for two lossless systems these equations conserve 
power, and for a lossy system (such as a nanotip) the effects of losses are treated exactly. By convention, the integrals 
appearing in the diagonal matrix elements are typically set to 1. 

For the waveguide and nanotip systems coupled over a length Lex: the exact single-photon generation efficiency will 
depend on the details of how the two systems are brought together and separated apart. In practice, for example, 
the two systems should be brought together slowly enough that the introduction of the waveguide does not cause 
significant back-scattering of the plasmon, yet quickly enough that this introduction length is small compared to 
the plasmon decay length. Furthermore, in reality the coupling region will not be a step of length Lex but will be 
characterized by some smooth transition. To avoid the many details associated with this introduction and separation 
and to approximately calculate the efficiency, we will consider an idealized system and make three assumptions; 

(i) The decay rates of the emitter are not affected by the presence of the nearby dielectric waveguide. In particular, 
the Purcell factors and error probabilities calculated earlier for the nanotip are unchanged. 

(ii) The radius of the nanotip is given by p{z) = ^/wz for z < zq and becomes constant, R = p{zo) — ^/wzq, for 
z>zo. For z>zo the plasmon mode solution becomes identical to the nanowire solution, and in particular has 
well-defined fcy which allows it to be easily mode-matched with the waveguide. It is assumed that coupling 
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between the nanotip and waveguide begins at z = zq, with the initial field amplitudes of the coupled system 
given by 

C^{zn) = (i-Pe{R)Y^ : (59) 
CgAzo) = 0, (60) 

where the effective plasmon excitation probability 1 — Pe{R) is already optimized for a given R over the nanotip 
curvature and emitter position. 

(iii) Eq. H56() exactly describes the coupling between the two systems in the region zq<z<zo + L^x- To estimate 
the probability of transfer from nanowire to waveguide after distance L^.^ when the two systems are once again 
separated, we project the total field of Eq. H55|l at z = zq + Lex into the waveguide mode. Specifically, the 
projected field amplitude in the waveguide in either of the degenerate modes i is given by 



C'proj,i 



(zo + Lex) = 2 J dp (ET(r) X H;,(r)) • z, (61) 



where the factor of 2 arises due to the normalization convention adopted here for the unperturbed modes, 
^Jdp (Eg,, X H;^, + X H,,,) = 1. 

Because of the symmetry, the projected field strengths |Cpi.oj,iP calculated above are equal for the two degenerate 
waveguide modes, and the quantity 2|Cpi.oj,iP then corresponds to the efficiency of single photon generation. Here the 
additional factor of 2 accounts for the mode degeneracy. This quantity takes completely into account the propagative 
losses of the plasmons, imperfect coupling between the nanotip and waveguide, and the Purcell factor of the nanotip. 

In Fig. we plot the efficiency of single photon generation as a function of i?, for both the nanowire and nanotip 
systems. For each R the plotted efficiencies have been optimized over all other possible parameters of the system. 
For the nanowire configuration, we have assumed that the resulting forward- and backward-propagating waves in the 
waveguide can be perfectly combined. In the figure we have also included points obtained by our BEM simulations of 
a nanotip. Here, we have taken the numerically optimized values of Pe and plugged them in as initial values for the 
coupled-mode theory above. It can be seen that the numerical simulations agree well with our theoretical predictions. 
We find that photon efficiencies of approximately 70% are possible for the nanowire, while efficiencies exceeding 95% 
are possible for the nanotip. In Fig. 03 we plot the optimal coupling length L^x, in units of Api, as a function of R for 
the nanotip {L^x for the nanowire should be twice that of the nanotip, to account for the transfer of the forward- and 
backward-propagating components of the emitted plasmon). It can be seen that the out-coupling to the waveguide 
can in principle occur quite rapidly, over length scales of a few Apj. 

The existence of an optimum R for photon generation can be intuitively understood. For smaller R the coupling 
between the emitter and plasmon modes can be quite large. However, these tightly-confined plasmon modes are 
accompanied by higher propagative losses and cannot be as efficiently coupled to the waveguide system. The coupling 
efficiency between plasmons and waveguide modes improves for larger R. For the nanowire, however, the larger radius 
results in weaker coupling between the plamson and emitter, while for the nanotip the accumulated propagative loss 
increases as the final radius R grows. 



VI. EFFECTS OF SURFACE ROUGHNESS 



In previous sections we have treated the problem of plasmon propagation on smooth nanowires and nanotips, 
taking into account inherent dissipative losses characterized by Im 62. In practice, however, these structures are not 
perfectly smooth, and the surface roughness can give rise to new scattering mechanisms for the plasmons. While 
the general solution for the fields in the presence of arbitrary roughness is a complicated problem, we calculate the 
effects in two limits. In Sec. IVI Al we calculate the losses on a nanowire due to radiative scattering in the limit of 
small roughness and zero heating (Ime2 = 0). Here the plasmons experience no inherent loss due to the metal but 
can receive momentum kicks from the roughness that cause them to scatter radiatively. In Sec. IVI 51 we calculate the 
effects of small roughness for a nanowire in the non-retarded limit, where radiative effects are ignored but the effects 
of increased dissipative losses are treated. While we explicitly treat only the nanowire case here, we note that the 
results obtained can also be incorporated into our model for nanotip losses via Eq. H54|) . 
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A. Radiative losses 



For simplicity we consider a wire with axial symmetry, but with a surface profile given by po{z) = R + p(^(z), where 
R is the average radius of the wire, C.{z) is some random function describing the roughness, and p is an expansion 
parameter that will be taken to equal 1 at the end. We will calculate in perturbation theory the radiated field scattered 
from the roughness, from an initial field corresponding to the fundamental plasmon mode for a perfectly smooth wire. 
Because of the symmetry, the only non-zero components of the fields remain Ep, E^, and 7?^, which will also have 
axial symmetry. As will be seen later, it suffices for now to consider only Ez, as the other components depend on 
Ez in a simple way through Maxwell's Equations. We proceed by breaking up the total field along z in region i into 
incident and scattered fields 

4°*"' = ^^+^^.' (62) 

where E'^ ^ is the z-component of the fundamental plasmon mode given by Eqs. (O and and further assume that 
the scattered field can be expanded in a power series 

oo 

i?L = Ep"4'- (63) 



In the following we will calculate the first-order scattered field E^ ^ . We make the ansatz that E^ I can be expanded 
in the form lisll 



/oo I 2 

d/^||i7o(/ii±p)^A(/i||)e"^"^ 

/oo I 2 

d/^ll Jo(W)-^i3(/i||)e^""^ (64) 
-oo "-2 

where each Fourier component is an outgoing solution of the wave equation with appropriate boundary conditions at 
p — ^ and p — oo, as derived in Eq. ljA4|l . From Eq. ljA4p one also sees that Ep, are determined completely once 
Ez is known. Using these relations, the total (incident plus scattered) fields E*"*^' and H*°*''^^ are straightforward but 
lengthy to write down, and are given to order p in Eq. (|F1I) in Appendix IfI 

The coefficients v4(/i||), B{h\\) are determined by enforcing continuity of the tangential fields at the boundary po{z). 
Specifically, we require that 

(1 TTitotalNl (7 TTitotalNl 1 ^ ' P dzP 

)\p=R+paz) - )\p=R+pazy^- 



rrtotall rrtotal I (P,K\ 

-"0,1 \p=R+pC{z) ~ ^4':^ \p=R+pC{z)' 

where t{z) is the unit tangent vector to the interface. These equations can be solved perturbatively by expanding 
them in p and solving at each order. It should be noted that the expansion should be done carefully, as dependence 
in p is contained not only in the fields given in Eq. IjFip but also in the surface profile po{z) = R + pC,{z) and tangent 
vector t. The 0{p'^) equation is trivially satisfied by the fundamental plasmon mode for a smooth nanowire. To solve 
to 0{p), it is useful to first introduce the Fourier transform of the surface roughness, 

/oo Jlj 
^e^h^aiH^). (66) 
-oo 

Using the Fourier transform Cih\\), the 0{p) equations become algebraic in Fourier space and have solutions (see 
Appendix IF|| 

m) = -^^-fZ-a^hl (67) 
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where fc|| denotes the unperturbed plasmon wavevector (in this section we take Im £2 = so that and C_i are purely 
real) . The scattering coefficients are complicated functions of and R and are given in Appendix IFl 

Physically, the equations above state that, to first order, the surface roughness contributes single momentum kicks to 
the unperturbed plasmon fields with a strength determined by the Fourier components of the roughness. From this 
point forward we set p = 1. 

We now consider some random surface profile such that 

(C(z)) = 

(C(z)C(^')) = <5^e~(^~^')^/'^^ (68) 

with corresponding correlations 

(C(fc)) = 

for the Fourier components. Physically 6 and a correspond respectively to the typical amplitude and length of a rough 
patch on the surface of the wire. It is also useful to define s — S/a as a typical "slope" to the roughness. To calculate 
the power radiated due to the surface roughness we will find the ensemble-averaged Poynting vector far from the wire. 
It is sufRcient to consider just the component of (S) oriented along p, given outside the wire by 



Sp = -^{E\T'HC°'-'), (70) 



where the fields £'5;°*'^' , H^°^^^ are given to first order by Eq. (|F1|) . The calculation of Sp simplifies further because the 
incident plasmon field decays exponentially away from the wire, and thus to lowest order only the first-order scattered 
fields will contribute to the Poynting vector at large p, which physically corresponds to the power radiated away to 
infinity. Specifically, the radiated power per unit area is given by 

Sp ^ -liEi'lHif) (p^^) (71) 



= ^ l°°^d^dh\'-^^^^^H^^^^^ (72) 

Substituting the solution for A{h\\) derived in Eq. (|67|l and using the correlations in Eq. I|f)9() . it is straightforward to 
evaluate the integral over and arrive at 

Sp - d;^||e-i-'(^|-fc|i)'/ii^ffo(/ii^p)i/^*(/^i^p) |/(/i||)|2 . (73) 



-fei 



In the expression above we have truncated the bounds of the integral to ±ki because we are interested in the Poynting 
vector far away from the wire, where only radiative fields j/iy | < ki should contribute. With knowledge of the Poynting 
vector it is then possible to find the dissipation rate of the plasmons due to radiative scattering, given by 

y. _ y 2TTpSp 

rad, rough ^ ^™ 1 r , ^ / / \ ni-.-,/ \i2 i-r-r/ m2' \'^) 

i/dpeo£ (e(p,c^)^) |E(p)| V MO |H(p)r 

The denominator on the right-hand side of the equation above can be identified with the plasmon energy per unit 
length. 

We first qualitatively discuss the behavior of ^J.g^(J ^.Q^gJl before deriving various limits more quantitatively. From 
Eq. (|73|l it is clear that F^ad, rough scales explicitly like 6^ or s^. Physically, this occurs because the lowest-order 
contribution to the Poynting vector far away from the wire is due to the combination of a first-order scattered electric 
field and first-order scattered magnetic field. In Fig. [T] the quantity rj.ad,rough/'S^'^ is evaluated numerically as a 
function of wire radius R and correlation length a/R, for a silver nanowire at Aq — 1 pm and ei = 2. We are 
particularly interested in the nanowire limit, when the plasmon wavevector fc||«C_i/i?. We see that for fixed R, 
the scattering reaches a peak for some particular value of a/R. More careful inspection reveals that the maximum 
occurs when aoci?/C_i oc Api. This result makes intuitive sense, since the characteristic momentum kick '^l/a that 
the plasmon wavevector fcy receives due to roughness must be on the order of C-i/R in order for the resulting 
wavevector to lie in the radiative range between —ki and fci. In the limit a/i?3>C_i, one observes an exponential 
suppression of scattering, due to the fact that the roughness has a very narrow momentum distribution and cannot 
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possibly contribute a large kick to . In fact, in this regime one physically expects for the plasmon wavevector to 
adiabatically vary with the changing wire radius. In the other limit a/i?<cC_i, the scattering also decreases, but 
with a polynomial dependence on R, as will be proven below. Here, the momentum distribution of the roughness 
becomes very wide, and thus the probability of receiving a kick that results in a final momentum between ±ki becomes 
quite small. Finally, for fixed slope s, it can be seen that the scattering decreases as R^Q at any correlation length 
a. This result is also easily understood, as the plasmon wavevector k\i becomes increasingly far-removed from the 
range of radiative wavevectors. In Table HI we calculate the scattering rates for wire sizes koR = 0.1,0.2,0.3 (or 
i?wl6,32,48 nm), for a few chosen roughness parameters. The scattering rates are given as a percentage increase in 
Im fc|| over the values for a smooth nanowire. It can be seen that strong suppression of radiative scattering occurs 
both for smaller R and when a is either much larger or much smaller than R, which confirms our earlier observations. 
Furthermore, it is evident that under reasonable parameters, the losses in the system are increased only slightly due 
to radiative scattering, around an amount of 10% or less. 

We now analyze more carefully the behavior of the radiative scattering in the nanowire regime. For concreteness, 
we will consider a field normalized by Eq. (|29|l . in which case the denominator of Eq. (|74(l becomes huj/AL. To simplify 
the expression further, we first note that since we are interested in the far field {p — > cxd), we can take the asymptotic 
limits of the Hankel functions in Eq. H73|) . Ho{hi±p)HQ* {hi±p)K — 2i/{7Thi±p). One can also derive an asymptotic 
relationship of as i?— >0 (see Appendix IF|) . which upon substitution yields 

rrad,rough « n'/'^-€cue,s'a' hl^e~^^'''W^h-^^0\ (i?^o) (75) 

<!> ^ ^°ff'\^ 4^(»g-i)-go(»g-i)- (76) 

From the equation above, it is clear that there are three distinct regimes of interest defined by the quantity 
a=k\\a = 2TTa/ Xpi^C-ia/ R, which characterizes the typical extent of a rough patch compared to the plasmon 
wavelength. In the limit a^l, one can approximate the exponential in the integrand of Eq. H75|) as a constant, 
which leads to straightforward evaluation of the integral, 

T. _4 5/2 2 Ao-R"\^ 3 f 

rrad.rough-gTl" / —^^^1 S (a<l) (77) 

Here, the noise spectrum of Eq. Ht)9|) becomes very wide and leads to an scaling of the dissipation rate. In the 
opposite limit q;3>1, the value of the exponential term becomes exponentially small, with a corresponding exponential 
suppression of the scattering rate. A more careful evaluation of the integrand yields 

Finally, one can show that for fixed, sub-wavelength R the radiative scattering is most significant when a^O{l). 
In this case, the exponential appearing in Eq. H75|) is neither exponentially small nor constant. However, one can 
make the rough approximation e~^^/^^° Kil — (l/4)a^(fc|| — /i||)^ to get an idea of the scaling in this regime. 

It is straightforward to show that the scattering rate has a maximum with respect to a at a«(12/5)^/^, with a 
corresponding maximum decay rate 

!r X l"?^!^ 5/2 2 f koRV , . 

max|ri.ad,rough| «— we/ s [-^j ■ (79) 

Again, the radiative scattering is most significant when the length scale a of the roughness is on the order of the 
plasmon wavelength, and the maximum scattering (for fixed s) decreases as i?— >0 due to the increasing mismatch 
between fcy and radiative wavevectors. 

We now consider the limits of validity of the derivations above, specifically considering the expansions made in 
Eq. ljF2p that are necessary for the perturbative method used here. The first of these expansions requires that 
\dC,/dz\<^l, which can be re-written as a condition on the slope, s<Cl. Physically this requirement states that the 
typical length of a rough patch be much larger than its typical height. The second line of Eq. (jF2|) requires that 
|A;ij_C|<Cl. In the nanowire regime this requirement is equivalent to (5<cApi, which states that the height of a rough 
patch must be much smaller than the plasmon wavelength. Finally, the third line of Eq. (|F2I) requires within 
the range of hi± that are appreciably scattered into. From Eqs. H()7|) and H()9|l . we see that the relevant range for the 
parallel component of the wavevector is given by fc|| — 1/a ^ ~ fc|| + 1/a, and thus the largest relevant transverse 
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wavevector is |/ii_L,max|~ inax_i<£)<i \^ei(ijj/cY — (k\\ + 9/a)'^\. In the nanowire regime, fc||«C_i/i?, there are two 
Umiting cases. The first is when the correlation length a is much larger than R, a^R, in which C&S6 |/i2_L max 
and |/ii_LC|<^l reduces to (5<CApi. In the other limiting case, R^a, one finds that |/ii_L,maxHl/'^ ^^'^ the corresponding 
requirement is given by s ^ 1. 

We finally note that while the radiative scattering goes like 6^ or , the relevant quantity for dissipativc (heating) 
losses due to roughness becomes S inside the wire. For this quantity the lowest-order correction to the smooth wire 
solution will come from a combination of a first-order and zeroth-order field. Thus one expects roughness-induced 
dissipative losses to contribute a decay term proportional to (5 or s, which for small roughness will dominate over 
radiative scattering. This correction will be treated in the next subsection. 



B. Non-radiative losses 



To study the effects of surface roughness on non-radiative losses, we will make one simplifying assumption and 
calculate these losses in the quasi-static limit. To do this we will proceed in a manner similar to that in Sec. IIII Al 
where we found the quasi-static fields associated with a smooth nanowire. Here the calculations for the fields yielded 
the presence of poles whose positions and widths give the real and imaginary parts of the wavevector fc|| . The case 
of a smooth nanowire was particularly easy to treat because of the translational symmetry of the system. A system 
containing surface roughness lacks such translational symmetry and therefore must be considered more carefully, 
but the calculation proceeds in much the same way. In particular, we will find expressions for the pseudopotentials 
$1 — <I>o + <&r and $2 that satisfy the necessary boundary conditions in the presence of surface roughness. We can once 
again find the positions and widths of the poles associated with the system, which are now altered by the roughness. 

We first write down appropriate expansions for $o.r,2(r, r')- The expansion for the incident component $0, given 
originally in Eq. Hll() . is slightly re- written here. 



$o(r, r') 



1 



47reoei |r - r'| 

^ (2- 5„,o) cos (m((/) -(/)')) / dh cos {h{z - z')) K,n{hp')I„^{hp) {p < p') 



m— — 00 



47r2eoei 

The functions Km{x), Im{x) are defined by 



. ^ e^W-*') / dhe'''^'-''^kM)L{hp) {p<p'). (80) 



KmJm{x) = i^™,/™(|x|). (81) 



We also break up $^,2 into Fourier components that satisfy Laplace's equation, and assume that these expressions 
hold up to the interface: 

1 roo 

$,(r,r') = — V e™(*-*') / dhe^'^^'-''^ Kr,{hp)a^{h), (82) 

$2(r,r') = V e-'(^-^') / dhe'^^'~''U^{hp)f3^{h). (83) 



m— — OQ 



To describe the surface roughness, we assume an interface with axial symmetry as before, Pq{z) = R + p(^{z), where 
the roughness profile ( satisfies the correlations given in Eqs. (|68|l and Ht)9|) . The coefficients a^, are determined 
by the boundary conditions, namely that $ and must be continuous at the interface: 



^i"-V*i(r,r')|,=«+,,(,) = ^2--V$2(r,r')|^^^^^^(,), n = ^£=£fc=. (84) 



Plugging in the expressions for $i given by Eqs. (|8()|l . H82(l . and H83(l . one can expand both boundary condition 
equations to 0{p^). Then, replacing C(z) in these equations with its Fourier transform given in Eq. one obtains 
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a set of coupled equations for arm Pmih) completely in Fourier space, which, unlike the case of a smooth nanowire, 
is not de-coupled in h. It is tedious but straightforward to show that, to 0{p^), this set of equations is given by the 
matrix integral equation 



Mh)+P I ^ah-q)Mh,q)+p' f ^ah-q-q')aq>2ih,q,q')+Oip^). (85) 



27r^' ^' >iy ■ ^ y (27r)2 

The matrices Mi and vectors are complicated expressions and are given in Appendix [HI We note, however, that 
in the case of no surface roughness (C = 0), the solution to the resulting equation Mo{h) ■ {am{h) f3m{h))'^ = vo(/i) 
reduces to that of a smooth nanowire. 

We now discuss how to solve Eq. H85|l in the presence of surface roughness, using the methods detailed in |43]. One 
might first consider expanding am,f3m in a power series of p, in a manner similar to the field expansion in Eq. H63|) 
for the case of radiative scattering, and then solving the ©(p""*"^) equations based on the ©(p") solutions. However, 
one expects that such a perturbative solution would simply yield poles for each higher-order correction with the same 
location as that of the unperturbed solutions , /3m' . Mathematically, this occurs because each calculation of the 
next correction involves an inversion Mfj~^(/i). On the other hand, physically we expect for the surface roughness 
to result in some shift of the pole that is not predicted by such a perturbative method We thus consider an 

alternate approach, in which we symbolically sum the perturbation scries in Eq. (|85|l to all orders and then only keep 
the lowest order result in p. Let us symbolically write Eq. (|85|l in the form 

{Mo + 5M)x = ia + Si, (86) 

where Mq and fo arc non-random matrices and vectors, respectively, SM is a random 2x2 matrix integral operator. 
Si is a random vector, and x is a column vector with components am, Pm- We now define the averaging operator 

Px={x), (87) 

and the operator Q = 1 — P. We can apply P, Q to Eq. H8t)|l to get 

PMo:ii + PSMyi = P(fo + (5f), (88) 
QMqx + QSMx - Q{io + Si), (89) 

which after some manipulation results in the set of equations 

{Mo + PSM){K)+PSMQx^P{{o + 6i), (90) 
Qx = (1 + Mo^QSMy^Mo^Qifo + Sf) - (1 + Mo^Q5Mr^Mo^QSM(x) . (91) 

One can then substitute Eq. (|91|l into Eq. 190|) and solve for (x) , in which case one obtains 

(Mo + ((1 + SMMo'Q)-^6M)) (x) = ((1 + 6MMo'Q)-\fo + Sf)). (92) 

We note that, unlike a perturbative expansion and solution for am,Prrn the equation above is thus far exact. Now, 
we assume that SM and Sf can be expanded in powers of p in the form 

SM = pSMi +p^SM2 H , 

Sf = pSfi + pHf2 + ■■■. (93) 

Comparing Eqs. H85|) and H93|) . we see that {SMi) = (^fi) = since (C) — 0. With this result, and utilizing the 
definition Q = 1 — P, one can proceed to expand Eq. (|92|l up to ©(p^), which yields (after setting p — \) 

{Mo + {SM2) - {SMiM^'SMi)) (x) = (fo) + {Sf2) ~ {SMiMo'Sfi). (94) 

Substituting the corresponding terms of Eq. (|85|l into the equation above and using the second-order correlations 
given by Eq. H69() . we find after simplifying that 

Mo{h) + dq (e-'^''^'''M2{K h, q) ~ e''^'^'^-'^^' ^'Mi{h, q)M-\q)A'h{q, h)) ^ ^ = 

Mh) + ^J dq(^e--''^'/^^2{h,h,q)~e-'''^'^-''^'/^Mi{h,q)M,\q)vi{q,h)) , (95) 
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where s = 5 /a. 

We now discuss the solution to aQ{h), which contains a pole corresponding to the fundamental plasmon mode 
TO = 0. When £2 is a negative real number, ao has a pole on the real /i-axis whose position gives the new, shifted 
plasmon wavevector fcy. When 62 has a non-zero imaginary component, ao will have a resonance feature along this 

axis whose peak corresponds Re fcy and whose width corresponds to Im fcy . A quick inspection of the equation above 
reveals that A:||i? = C-i{ei, s,a/ R) is a constant that depends only on the quantities ei,s, and a/R. Unfortunately, 
because of the complexity of Eq. H95() it is difficult to derive other scaling results for even in limiting cases. However, 
Eq. H95|l can be solved numerically. In practice, for known parameters, the matrices M and vectors v can be readily 
evaluated over some range of h, from which the solutions to the system aojPoih) in that range immediately follow. 
The resulting resonance in ImQ;o(ft.) as a function of h is then fitted to a Lorentzian, with its peak giving the shifted 
wavevector Refcy and its half-width giving ImA:||. In Table ITll we give the resulting losses and wavevector shifts for 
a few roughness parameters, as calculated through Eq. H95|l . Again the numbers that we have used are for a silver 
nanowire at Aq = 1 /im and ei = 2. The shifts in ReA;|| (or equivalently, ReC_i) and increases in the loss parameter 
Im fc||/Re /c|| (or Im C_i/Re C_i) are given in terms of the percentage increase over their values for a smooth nanowire. 
Again it can be seen that for reasonable parameters, surface roughness adds only a moderate amount of loss to the 
system. 



VII. CONCLUSIONS AND OUTLOOK 



We have demonstrated that the subwavelength confinement of guided plasmon modes on conducting nano-structures 
leads to strong coupling between these modes and nearby emitters in the optical domain. This strong coupling leads 
to large effective Purcell factors for emission into the plasmon modes, which are limited only by heating losses in the 
conductor. While losses prevent the plasmon modes from being useful photonic carriers of information, we have shown 
that they can be efficiently out-coupled, e.g., to a dielectric waveguide. We estimate that single photon generation 
efficiencies exceeding 95% are possible in such a tiered system. Finally we have analyzed the effects of plasmon 
scattering due to moderate surface roughness on these nano-structures. 

Rapid advances in recent years in fabrication techniques for nanowires ^5^, '5l] , nanotips |52| , and sub- wavelength 
dielectric waveguides ,,53„ .54i| puts such a system in experimental reach. Quantum dots or single color centers might 
serve as physical realizations of solid-state emitters, which could be used to achieve strong-coupling cavity QED and 
quantum information devices on a chip at optical frequencies. It is also interesting to consider real, individual atoms 
interacting with nanowires and the challenges associated with constructing nanoscale traps. These traps might in 
part be formed by the plasmon fields themselves. 

We emphasize that the physical mechanisms that lead to strong coupling are not restricted to the nano-structures 
considered here but can be quite a general feature of the plasmon modes associated with sub-wavelength conducting 
devices. It is thus likely that the efficiencies calculated here are not fundamentally limited but can be further improved 
by proper design. Photonic crystal-like structures for plasmons js^, for example, may be a promising approach to 
achieve tight confinement while simultaneously reducing losses. Similar schemes may also help to improve coupling 
between the plasmons and dielectric waveguide modes. Such approaches are likely to improve the performance of 
plasmon cavity QED even further. 

The authors thank Atac Imamoglu for useful discussions. This work was supported by the ARO-MURI, ARDA, 
NSF, the Sloan and Packard Foundations, and by the Danish Natural Science Research Council. 



APPENDIX A: GENERAL THEORY OF ELECTROMAGNETIC MODES OF A CYLINDER 

The solution to the electromagnetic modes of a cylinder has been known for quite some time |2^ (2^ and is briefly 
derived here. 

We consider a cylinder of radius R of dimensionless electric permittivity £2, centered along the z-axis and surrounded 
by a second dielectric medium £1. For non-magnetic media the electric and magnetic fields in frequency space satisfy 
the wave equation 

v'{|3}4«"){hS}-- <-) 

The solutions to Eq. I|A1|I are perhaps most easily derived by first finding scalar solutions of the equation and 
then constructing vector solutions. Working in cylindrical coordinates, scalar solutions of Eq. HAljl satisfying the 
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necessary boundary conditions take the form ■0ioci?m {ki±p) e™'^+*'^ii^ and V'20c Jm {k2±p) e*'"'^+*'^ii^ outside and inside 
the cyhnder, respectively. Here J,n and are Bessel functions and Hankel functions of the first kind, respectively, 

while ki± = ^kf — /c^ and ki — uj^/ei/c. J„i is well-behaved at p = 0, while Hm{x)'^e" for large x satisfies outgoing- 
wave conditions. It is easy to verify that two independent vector solutions to Eq. IjAljl are given by v.^ = t-V x 



and Wi = -^V x Vj. The curl relations of Maxwell's Equations then imply that E and H must take the form 



E.(r) 
H.(r) 



6-iW.,(r), 
-h (aiWi(r) -) 



aiVj(r) 
i 

LOpQ 

where ai,bi are constant coefficients. Expanding out these expressions in detail 

ik\\kij_ 
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H.(r) 
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ki± I iTT-kw 
h ' klp 



{h±p) (j) 



k? 



aiF-j^{ki±p) + - — b,F,^„^{ki^p) p 



hp 



mk\\ 



h± 

aiF^^rn{hi_p) + -^biFl ^{h_Lp) 



-[^aiFt^rnih±p)z 



(A4) 



where = Hmix) and i^2,m(a;) = Jm{x). 

Up to this point Ui, hi are arbitrary coefficients, whose relationship becomes fixed by imposing boundary conditions 
between the cylinder and surrounding dielectric. Requiring that the tangential components Eip, E^, H^, of the 
fields be continuous at the boundary results in a linear system of four equations, which we write in abbreviated 
matrix form as M{ai 02 61 62)"^ = [s^. A non-trivial solution for the fields requires that detM = 0, which after 
some work simplifies to the mode equation given in Eq. 

One special case of interest is that of a TM mode with no winding (m = 0). The component of H along z by 
definition vanishes, which implies that the coefficients in Eq. l|A4p vanish. The condition det M = is then 
significantly easier to evaluate in this situation. In particular, ai — implies that the field components E^ and 
vanish, and continuity of the remaining tangential field components Ez and at the boundary requires that 



Hoiki^R) 
kii_H'^{kiA^R) 



Ulfj.0 



'^Mk2^R) 

k2±4ik2±R) 




(A5) 



Setting the determinant of the above matrix equal to zero immediately yields the mode equation of Eq. |5J , and it is 
also immediately seen that the ratio of the coefficients 61,2 must be given by Eq. 



APPENDIX B: DERIVATION OF CUTOFF FOR HIGHER-ORDER MODES 



In this section we show that to a very good approximation, a nanowire essentially supports a single, fundamental 
TO = plasmon mode. In particular, for all higher-order plasmon modes |m|>2 a cutoff wire size i?cutofr exists below 
which such modes cannot exist, while the |m| = 1 plasmon modes exhibit an exponential growth in their mode volumes 
as i?— s-0. For simplicity, we will assume in this section that we are dealing with a lossless system (Ime2 = 0). 



1. Behavior of |mj>2 modes 



We are interested here in the behavior of the |to|>2 modes near cutoff, which is characterized by a small deviation of 
the plasmon wavevector fcy from y/eiui/c (see Fig.EJ- To simplify algebra in the derivation of -Rcutoff: from this point 
forward we make the mode equation ^ dimensionless by setting uj/c = 1, and we will assume that m is positive (the 
case where to is negative follows this derivation with a few minor modifications). Furthermore, it is useful to define a 
small quantity S = k\\ — where we specifically consider the positive fcy solution. On physical grounds, any mode 
with positive fcy must satisfy 5>0, because if k^^ < the fields outside of the wire would be radiative in nature and 
implies that the system is continually radiating energy out to infinity without a source. It follows that any value of R 
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where 6 = becomes a solution to Eq. for some m then corresponds to a critical point in behavior, and specifically 
is a cutoff beyond which modes cease to exist for that given m. To find this R = ^cutoff j it is useful to expand the 
two sides of Eq. Q in 6. We will find that both sides have contributions to these expansions that are divergent at 
(5 = (terms that behave like S~", where n > 0), and we will show that, for m>2, there exists one value of R that 
equates these two divergent contributions; i.e., S = satisfies the mode equation at this particular value R = i?cutoff- 
It is straightforward to show that the divergent contribution to the expansion of the left-hand side of Eq. |^ is 
given by 

LHS - - ^^^^^=^---^ + 0{S ). (Bl) 

To expand the right-hand side, we first note that the quantity {^1 k2A_)Jl^{k2^R) I Jm{k2^R) = 
{l/y/e2-~ei)J,n{y/e2~~eiR) + 0{S^) is well-behaved near 6 = 0. Here we have defined Jm{x) = J'^{x) / Jm{x) . Then, 
using the identity 

Hrniix) = —Lri^M, (B2) 
where K„i{x) is a modified Bessel function of the second kind, and the expansions 

Kraix) = ^^^^ ( - ) - '-^^—^ ( - ) + 0(a:^-™) (m>2), (B3) 




k,^ = i \ \I26^, + -^^^ + 0(<55/2)J (B4) 

= iKii, (B5) 

it is tedious but straightforward to expand the expression {l/ki±){Hj^{ki±R)/ Hm(ki±R)) as well. Performing these 
expansions and simplifying, one finds that 



''''^-m-AR^6^W^5^ 2R^e^{e,~e2)6 ^^^^^^ ^""'^ 

Comparing Eqs. ljBl|) and ljB6|l . we see that 5 = is a solution provided that these terms are equal to 0{6~^), i.e., 



TO ei + £2 _ Rei ^ ^(ei + £2)>/m(Ve2 - £ii?) 
i? £2 — ei m - 1 Vei - £2 

The solution R = i?cutoff to Eq. (|ij7|l gives the cutoff wire size below which the mode m cannot exist. In the regime 
of interest (ei > 0,62 < 0,ei -1-62 < 0), the first and second terms are positive while the third term is a negative 
function (for R > 0) that behaves like —1/R for small R and approaches a constant for large R. It can be seen then 
that a solution exists for any m>2, which establishes that these modes are indeed cut off in the nanowire limit. 



2. Behavior of m = 1 mode 



For simplicity we will assume that to = 1, as the case of to = —1 follows this derivation closely. The case of m = 1 
must be studied separately because the expansion of Km{x) given in Eq. (|B3|I only holds for to>2. The different 
asymptotic scaling of Ki{x) leads to unique behavior of the m = 1 mode in the nanowire limit. In particular we will 
show that this mode does not strictly have a cutoff size, but that fcy — > ^Jel exponentially in the limit -R^O. In turn, 
the magnitude of fci^ becomes exponentially small, which corresponds to an exponential growth in the spatial extent 
or mode volume. 

Again defining 6 = k\\— ^ftl, we are interested in finding an approximate solution to Eq. in the limit of small 
6. We proceed by expanding both sides of the equation as a series in the small parameter. The expression for the 
left-hand side given by Eq. IjBljl remains valid for to = 1. For the right-hand side, we anticipate that both the 
quantities ki±R and ki± will be small as R—^0 (these assumptions can be checked for consistency at the end of the 
calculation), and we thus expand around ki±R = the term 

1 H[{kij_R) _ 1 H'Sn.^R) ^gg^ 



fci^ Hi{ki±R) iKi± Hi{iKi±R) 



^ ^ ' -^(7 + log^)K,^R + 0{Ki^R')], (B9) 



iKl± \Kl±R 
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where 7Ri0.577 is Euler's constant. Here we have used Eq. ljB2p to convert Hm{ix) to Km{x) and the expansion 



X \2 2 4 2 _ 
We now assume that ki±R is small enough that logKij^i?|, such that 



~ ' - zKi^fllog (Bll) 



fcix Hi{ki±R) iKi± \Kij^R 

Furthermore, having assumed that ki±^ (and by extension, kij^) is a small quantity, we can now expand the expression 
above in terms of 5 using Eqs. HB4() and (|ij5|l . Making this substitution, and after a bit of algebra, one finds that the 
expansion of the right-hand side of Eq. is given by 

1 e, + ie,-2R\,[e^-e,)\og{5R■'^,/2) 
^^^-4i?5^ + AR^^,[e,-e,)5 " ^^^'^ 

Finally, equating the left- and right-hand sides to 0{5^^) gives the solution 

exp —\. (B13) 



i?2^ "V i?2ei(e2-ei 
It follows that in the nanowire limit, 

Ki± - (kl-ei^^ (B14) 
« (2(5V^)^/' (B15) 

Eqs. IjBlSp and (jB16|) indicate that the ra — \ plasmon mode does not have a cutoff in the nanowire limit, but 
instead that its longitudinal wavevector approaches y/ei exponentially, with a corresponding exponential increase in 
its transverse extent (~1/kij^) and mode volume. It is therefore well-justified to say that this mode is effectively cut 
off, as the coupling strength to this mode becomes strongly suppressed as i?^0. 



APPENDIX C: RADIATIVE AND NON-RADIATIVE DECAY RATES NEAR A NANOTIP 

Here we derive more carefully the expressions given in Eqs. (|47|l and (|48|) for the radiative and non-radiative 
spontaneous emission rates near a nanotip. 

To calculate the radiative rate, we should consider our expression for in Eq. H42(l in the far-field (large v) limit, 
where the Km{qv) terms in $r decay exponentially with v. Because of this exponential dependence at large v, to good 
approximation it suffices to expand the terms am{<l), Jm{Qu'), Km{qv') around q = Q. The only non-trivial expansions 
occur for the terms am{q) and are given by 

«o(9) = \ (l-'Aq\l + 0{q% (CI) 



2 V ei 

= \'-^q'vl + 0{q% ... (C2) 

Z ei + £2 

These expansions allow for exact evaluations of the integral. It can be verified that the dipole contributions to 
originate from the m = 0, 1 terms in the sum, which are readily found to be 

1>^^°Hr,r') « -r^^^.{l-'A%^,yn^ + 5m, (C3) 
e-Hr,r') . r^cos(0-0')^^^^^. (C4) 
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Here 5$ is a complicated function, but most importantly contains no dependence on r'. Recalling that the pseudopo- 
tentials derived above correspond to a point charge source, we can immediately obtain the potentials due to a dipole 
Poe^*"* at r' by applying the operator (po • V) to these expressions. In parabolic coordinates the gradient operator 
is given by 

, ^ J n—+v— +— 0— , C5 
and for a dipole located on the z-axis {u' = 0), wc find that 

(m=l) _ 1 ei - £2 vl (po - «(P0 ■V))-V , . 

dip,r ~ 47reoei ei + ea ^3 ■ ^ > 

From these expressions one can immediately identify the induced dipole moments in the nanotip, 

5y> = -^^'o^(l-l)' (Pop) (C8) 

<5p = ^iPo^^4' (Po^^) (C9) 
ei + £2 w 

and arrive at the radiative decay rates given in Eq. H47|) . 

The leading term for the non-radiative decay rate is found by calculating the divergence in the reflected field Er(r', r') 
as v^vq. The reflected field = — V(po • V')$r is in general difficult to evaluate, but simplifies considerably for a 
dipole located on-axis {u' — 0) due to the presence of the Jm{qu') term in <i>r, given in Eq. The operation V 

causes terms like Jm(0) and J^i^) to appear in E,., which are non-zero only when m = and m — 1, respectively. 
This immediately leads to the expressions 

„2 roc 3 

Po-E,(r',r') = --^ / dq ^a^{q)K!{qv'), {po ± z) 

Po-E,(r',r') = / dq ^ao{q)KUqv'), (po || 5) (CIO) 

which were given in Eq. (|46|l . Examining further the solutions to ao,i, it can easily be shown that their asymptotic 
expansions in the limit qvo ^ 1 take the form 

ao,i(q)«-^i^e2«'"'. ((j«o»l) (Cll) 

TT ei + 62 

At the same time, in the limit qv' ^ 1 the behavior of is given by Kf{qv') w {tt /2qv')e^'^'^'" , and thus as w'— >wo 
the integrands of Eq. (jClOjl exhibit very long tails due to the presence of terms ^e^^"^'" at large q. The tail is 
the origin of the divergence that we expected on physical grounds. Using these expansions as well as the fact that 
the decay rate is proportional to Im (po • E), the integrals can be evaluated exactly and yield the non-radiative decay 
rates given in Eq. (gSJ. 



APPENDIX D: BOUNDARY ELEMENT METHOD 

Our numerical implementation of the boundary element method (BEM) closely follows the method derived in |4lj | . 
Here we briefly outline the main ideas of BEM while referring the reader to |41,] for more details, and we discuss the 
key elements of our implementation. 

We assume that our system contains a set of known, time-harmonic source charges and currents Pextjjext in the 
presence of some scattering dielectric body whose surface is denoted S (although we discuss one body here, BEM is 
easily generalizable to treat multiple scatterers). In the case of interest, S represents the surface of a metallic nanotip, 
while the external source corresponds to an oscillating point dipole poe"*'^* at some location r'. For simplicity we 
also assume that we are working with non-magnetic media, and we denote by (j ~ 1, 2) the dimensionless electric 
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permittivities outside and inside S, respectively. The underlying principle behind BEM is that the scalar and vector 
potentials (f)j (r) and Aj (r) in each region can be written (in the Lorenz gauge) in the form 

0,(r) = / dr'G,(r-r')pcxt(r') + ^^ / dsG,{r ~ s)a,{s), (Dl) 

A,(r) - g 1 dr'G,(r-r')jext(r') + £^dsG,(r-s)h,(s), (D2) 

pikjr 

Gj{r) = (D3) 

where Gj is the Green's function in a medium of uniform ej, and kj = y/ej{uj/c). Physically, the equations above state 
that the fields in region j can be described as a result of the combination of the external sources and some effective 
surface charge and current distributions cr^, hj on S. In general, these effective distributions do not have physical 
significance; for example, they do not correspond to actual charges and currents, and the distributions in region 1 
and region 2 are not necessarily equal {e.g., (Ji{s) ^ (72 (s)). The values of crj,hj are not known initially, but a set 
of linear integral equations for these distributions results from enforcing various boundary conditions for the scalar 
and vector potentials at S. In particular, 0, A, D^, and H|| must be continuous at the boundary. To calculate the 
distributions numerically, if the boundary 5' is finite, one can mesh up the surface into a finite number of grid points. 
Assuming that CTj, hj are constant over each grid point, the linear integral equations become a set of linear equations 
in the values of aj , hj that can be solved straightforwardly. Once these distributions are known, the potentials and 
then the fields E, H can be calculated. 

In our problem of interest, we assume that the dipole is located on the z-axis and oriented along z, while the 
nanotip is described by a paraboloid of revolution around the z-axis. Due to the axial symmetry of the system, BEM 
simulations are advantageous because one only needs to calculate the unknown distributions along one dimension 
instead of over the entire two-dimensional surface S. At the same time, the source is a dipole oscillating at constant 
frequency, and thus the external charges and currents are calculated quite easily. In BEM (at least in the current 
formulation), it is necessary that the nanotip surface S be finite, and we implement this numerically by tapering and 
rounding off the nanotip far from the region of interest. In general, any termination can result in some back-reflection 
of the guided plasmon, and this results in some small oscillations of the fields due to interference with the forward- 
propagating plasmon, as barely seen, e.g., in Fig. 01 In our simulations, the reflected amplitude is kept to within a 
few percent. Very flne meshes were used to ensure accuracy; in most of our simulations, for example, the spacing 
between points in the regions of constant R was approximately Api/400. 



APPENDIX E: DERIVATION OF COUPLED-MODE EQUATIONS 

In this section we derive the equations of evolution for two electromagnetically coupled systems based on Lorentz 
reciprocity. 

First we derive the Lorentz reciprocity equation generally. Assuming non-magnetic media, suppose that 
{E]^(r), H]^(r), e]^(r)} and {E2(r), H2(r), e2(r)} separately satisfy Maxwell's Equations. At this point the systems 
1, 2 and their field solutions are not necessarily related to each other at all. In the following we assume that all fields 
E(r, i) = E(r)e^*'^*,H(r, t) — H(r)e~*"* have harmonic time dependence. Using the vector identity 

V- (ax b) =b- (V X a) -a- (V X b), (El) 

and the curl relations of Maxwell's Equations we can write 

V • (El X H^) = H^ (V X El) - El • (V X H^) 

= • (iWA^oHi) - El • (ic^eoe^E^), (E2) 

and similarly 

V • (E* X Hi) = Hi ■ (V X E*) - E* • (V X Hi) 

= Hi • {-iuj^ioUD - • (-itJeoeiEi). (E3) 

Adding up Eqs. ljE2|l and (jE3|l yields the equation for Lorentz reciprocity, 

V • (El X H^ + E5 X Hi) = zweoEi • E^(ei(r) - e^r)). (E4) 
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We now derive coupled-mode equations for two waveguides based on the Lorentz reciprocity equation above. This 
derivation closely follows that of 1431 • We emphasize that the nature of the waveguides can be quite general, e.g., they 
can be any type of normal dielectric or plasmon waveguide. We let the indices fi,iy = a,b refer to the system consisting 
of waveguide a without the presence of system b, and b the system consisting of waveguide b without the presence of 
a. We also assume that the surrounding dielectrics for systems a,b are the same, i.e., e^(r — oo) = ej(r = cxd), and 
that the waveguides are co-propagating along the z-direction. It is assumed that the total electric field for the system 
consisting of waveguides a and b together can be written as 

E^(r) = CAz)Kir), (E5) 

iy—a,b 

with a similar expression for H. That is, we assume that the total field can be written as a linear superposition of 
the unperturbed modes of systems a, b. For the case where systems a, b each have one allowed mode, the index v 
refers to these unperturbed modes. In general, when a,b have allowed unperturbed modes, is understood to 
be an index that covers all of these modes. We can derive exact equations of motion for C^{z) by using Eq. (|E4p . 
Specifically, we will let the index 1 = T in Eq. (|E4|I refer to the total fields Ey(r), firp{r) and the dielectric profile of 
the combined system eT(r), while we will let the index 2 = /i refer to any one of the allowed, unperturbed modes of 
systems a, b. Substituting this into Eq. I|E4(I yields 



V • (E^ X H* + E* X H^) = tueo^^ ■ E* (e^(r) - e* (r)), (E6) 



u—ajb 

Applying Stokes' Theorem to this result gives 



§-,T. J d'p{a{z)E,xH; + C,{z)E;xIl,)-z = zujeo ^ / E, • E^ (e^(r) - e;(r)) . (E8) 



]y—a,b u—a.h 

The left-hand side can be further simplified, 



= a; E / rfV(a(2)E, xH*+a(2)E;xHj.z (E9) 

y—a.h 

E ^ (/ ^'z' X h; + e;i X Hj . 2\ + c.{z) ("^ / d?p (e, x h; + e* x hj • z\ (eio) 



E ^ (/ '^^P ^^'^ X H; + E* X HJ . + [lue^ J <f p E, • E;(e,(r) - e;(r))j , (Ell) 



where we have applied Stokes' Theorem on Eq. ljE4|) to get the last line. 

Substituting Eq. (|Elip back into Eq. ljE8|) yields a set of Na coupled, first-order differential equations: 



^^P.m(^) = -ioJ^oY.^.{z)K,,{z), (E12) 

P.M ^ I dV(E,xH;+E;xHJ.f, (E13) 
K,^{z) = J dVE,-E;(e,(r)-ej,(r)). (E14) 
We emphasize that these coupled- mode equations are exact within the ansatz given by Eq. HE5() . 

APPENDIX F: RADIATIVE SCATTERING 



Suppose that in the presence of roughness, the first-order scattered field is given by Eq. (|64|l . Using the 
expressions derived in Eq. ijJJ) for the unperturbed, incident plasmon field, and letting denote the unperturbed 
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plasmon wavevector, the total (incident plus scattered) fields to first order in p are given by 
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(Fl) 



The boundary condition equations in Eq. H65(l can be solved by plugging in the fields above, carefully expanding 
the equations as a power series in p, and then solving for each order of p, utilizing the expansions 



i = z+p^p + Oip^), 

F^,m{h±Po) - F,^rnih±R)+pCh^Fl^ih^R) + Oip^), 
F^^^ (h.^po) = F,^^{h,^R) + pCh,^Fl„^{h,^R) + 0{p^), 



(F2) 



where Fi^^nix) — iJm(x) and F2^m{x) = Jm{x). The resulting 0(jP) equations are trivially satisfied by the plasmon 
fields of a smooth nanowire, while the 0{p) equations are found to be 
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Here we assume that the metal inherently has no losses, i.e., Ime2 = 0, such that fcy is purely real. Then, by plugging 
in the Fourier transform of C,{z) given in Eq. (|66|l . the equations above become purely algebraic. It is tedious but 
straightforward to show that the solutions are given by Eq. H67II . with the coefficients f{h\\) and g{h\\) defined via 



/(/ill) 
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biH'^iki^R) - 
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(F4) 



To evaluate the expression for Fj-ad rough Eq. H74() . it is convenient to normalize the plasmon fields using Eq. H29|l . 
in which case the denominator of H74|l becomes Tilo/AL and the normalization coefficient for the field outside the wire 

is given by bi~. — - ° ^ ' , as derived in Sec. IIII Bl Then, using the relationships knK.C / R, ki±fiiiC-i/ R one can 



y eoVC'tiL' 

calculate the leading terms of as i?— >0, 



fih) 



h2i_bi{C^ilRf^ 
—2ikih2± /nhi±R' 



(F5) 
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where we have defined (j)={b2/bi)jQ (iC-i) — -ffp (iC_i). In the equation above we have explicitly given the leading 
terms of the numerator and denominator of f(h\\). The ratio 61/62 is given in Eq. Q), which in the nanowirc limit 
results in the simplification of (j) given in Eq. 1)76(1 . 



APPENDIX G: NON-RADIATIVE SCATTERING 



The elements of the matrices Mi and vectors appearing in the matrix integral equation (|85|1 in the presence of 
surface roughness are given by 

A/o(M = r\^''.?mV (Gl) 



heiK'„^{hR) ~he2l'ra{hR) 

M (h = ( lK'm.{<lR) -QlLilR) ] pi(h-q)z' /pr,N 



2 - , 



M2ih,q,q') = [ Mr^'^^l ^IfF'^^J^'^'"'^''' 
V M^\h, q, q') M^^h, q, q') J 

Ml\h, q, q') = ei (^kZiqR) + q\h - q - q')k'^{qR) + \qq' {h - q - q')k'^{qR)^ , (G4) 
Mf{h, q, q') = -62 (trCi^R) + - q - q')iUqR) + ^QQ' (h - q - q'KMR)) , (G5) 



Mh) = - ( Ij:^ ] KM), (G6) 



-('^-^ ^ - ( er^iqR)lit,)u,R) ^ ^-^^y'''-''^ («^) 



V2(/i,g,g') = - ( 3. , \J'n{<lR) ^ W„(gp')e'^''~«^"'. (G8) 

\\r::[qR) + q\h-q-q')r^{qR) + hq'{h-q-q')I'^[qR) 
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FIG. 1: a. An emitter coupled to a nanowire is optically excited and decays with high probability into the plasmon modes of 
the nanowire. A single photon source is created by evanescently coupling the nanowire to a nearby dielectric waveguide over 
a length Lex- The single photon source can potentially be uni-directional, e.g., by capping one end of the waveguide with a 
reflective surface, b. An internal-level scheme that allows for shaping of the outgoing single photon pulses. An emitter that 
starts in state |s) is coupled to excited state |e) via a time-dependent external field Q.(t). We assume that the excited state 
|e) is coupled to state \g) via the plasmon modes, causing |e) to decay into \g) with high probability, while simultaneously 
generating a single photon in the plasmon modes. The shape of the photon wavepacket is determined by Q,{t). c. A similar 
scheme for single photon generation using an emitter coupled to a nanotip instead of a nanowire. Note that this scheme is 
naturally uni-directional, as the generated plasmons propagate in a single direction. 




FIG. 2: Allowed plasmon modes fcy as a function of R for a silver nanowire embedded in a surrounding dielectric ei = 2, for 
frequency corresponding to a vacuum wavelength Ao = 1 ^m and room temperature. The fundamental (m = 0) mode, in 
black, exhibits a 1/R dependence, while all other modes are effectively cut off as i?— >0. Inset: the propagative losses for the 
fundamental mode, characterized by the ratio Re fcy/Im fcy, for the same parameters. 



Roughness parameters 


fco-R = 0.1 (7?«16 nm) 


kaR = 0.2 (i?«32 nm) 


koR = 0.3 (i?«48 nm) 


a = O.liJ, 5 = 0.05i? (s = 0.5) 


0.09% 


0.5% 


1.4% 


a = O.li?, 5 = Q.IR (s = 1) 


0.4% 


1.9% 


5.6% 


a^R,5 = 0.05J? (s = 0.05) 


0.9% 


4.5% 


12% 


a = 5i?, (5 = 0.05i? (s = 0.01) 


2.8% 


8.0% 


10% 


a = im, 5 = Q.IR (s = 0.01) 


7.0% 


14% 


16% 


a = 2m, 5 = O.IR (s = 0.005) 


0.9% 


2.9% 


3.8% 


a = 257?, 5 = Q.IR (s = 0.004) 


0.3% 


1.3% 


1.8% 



TABLE L Losses due to radiative scattering off of surface roughness for nanowires of varying sizes and roughness parameters. 
The scattering rates are given in terms of the percentage increase in Im fc|| that one would expect over the values for a smooth 
nanowire. 
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FIG. 3: a. Solid line: Probability of error, Pe = 1 — ^pl/(r' + Tpj), in which an emitter fails to emit into the fundamental 
plasmon mode for a nanowire, plotted as a function of R and optimized over the emitter position. Dashed line: optimized Pe 
vs. curvature parameter w for a nanotip. Dotted line: effective probability of error, Pe = 1 — rpi(-R)/(r' + Fpj) for emission 

into a nanotip and successful propagation to final radius R. Solid points: effective error probability Pe for a nanotip, calculated 
numerically through boundary element method. Inset: same plot, zoomed in near RjW — 0. h. Contour plot of log^Q Pe for a 
nanowire, as functions of R and d/R. c. Contour plot of logj^g Pe for a nanotip, as functions of w and d/w. 



Roughness parameters 
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270% 



TABLE II: Losses and wavevector shifts due to non-radiative scattering off of surface roughness for nanowires with varying 
roughness parameters. The shifts in ReC_i and changes in loss parameters ImC-i/ReC_i are given in terms of percentage 
increase over the corresponding values for a smooth nanowire. 
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a) |Re(ExH*)| b) |Re (ExH*)|/r^^^^| 

k„d=o.oo2, r,^,„/r„=3.4xi 0^ k^d=o.oo2, r,^,^,/r^=3.4xi 0^ 




FIG. 4: Numerically calculated fields due to a dipole emitter near a conducting nanotip, obtained by boundary element method, 
a. The energy fiux |Re(E x H*) | , in arbitrary units. The position of the emitter is denoted by the blue circles, while the boundary 
of the nanotip is given by the dotted lines . The plots shown are for a final nanotip radius of koR = 0.3, curvature parameter 
kow = 0.022, and emitter positions kod = 0.002, 0.2, 0.7. It can be seen that both the total spontaneous emission rate F^otal ^^'^ 
the emission rate into plasmons increase as the emitter is brought closer to the nanotip. b. The quantity |Re(E x H*)|/r^-Q^-g^[, 
for the same parameters. This quantity is proportional to the energy fiux normalized by the total power output of the emitter. 
The kod = 0.002 plot is mostly dark, indicating that most of the decay is into non-radiative channels. The kod = 0.2 case is 
characterized by bright spots along the entire edge of the nanotip, which indicates efficient plasmon excitation. The kod = 0.7 
case exhibits the typical lobe pattern associated with radiative decay of a dipole. 
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FIG. 5: Wavevector fcy of the fundamental guided modes of a cylindrical dielectric waveguide with core permittivity ec = 13 
and surrounding permittivity ei = 2, plotted as a function of core radius Rg. 
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FIG. 6; a. Optimized eSiciencies of single photon generation vs. R. We have assumed that couphng to waveguide modes other 
than the fundamental mode is negligible, i.e., the waveguide is effectively in the single-mode regime. Solid line: theoretical 
efficiency using a nanowire. Dotted line: theoretical efficiency using a nanotip. Solid points: nanotip efficiency based on 
boundary element method simulations, combined with coupled-mode equations, b. Optimal coupling length Lex for a nanotip 
as a function of R. Here Lex is given in units of the plasmon wavelength Apj at that particular R. 




FIG. 7: The plasmon dissipation rate due to radiative scattering off of surface roughness, F^^^ rough/'^^'^i functions of wire 
radius R and correlation length a/R. The numbers are calculated for a silver nanowire at Ao = 1 fJ-m and ei — 2. 



